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The  mechanism  of  interaction  between  a  blast  wave  and  a  diffu¬ 
sion  flame  is  studied  with  the  objective  of  correlating  data  obtained 
in  shocktube  tests.  Based  on  a  scrutiny  of  the  interplay  between 
physical  transport  processes  and  chemical  kinetic  reaction  processes, 
a  Damkohler  number  P  is  identified  to  represent  the  inverse  of  the 
blast  strength.  Also  identified  is  a  fire  strength  parameter  q  which 
defines  the  fuel  character.  On  a  map  of  P(q),  regimes  in  which 
extinction  is  expected  are  distinguished  from  those  in  which  the  fire 
will  sustain  the  impact  of  the  blast.  Where  charring  fuels  are  in¬ 
volved,  the  role  of  glowing  combustion  in  the  reflashing  of  an  extin¬ 
guished  flame  is  examined. 

The  hypothesis  is  tested  with  the  experimental  data  obtained  by 
Backovsky,  Martin  and  McKee  at  SRI  International  to  reach  conclusions: 
(a)  liquid  fuel  (i.e..  Class  B)  fires  are  consistent  with  the  hypo¬ 
thesis;  (b)  barricades  placed  upstream  of  the  fire  tend  to  produce 
recirculation  of  flow  and  as  a  result  increase  the  flow  residence 
time  to  enhance  the  stability  of  the  fire;  (c)  the  barricade  heights 
for  which  data  are  now  available  are  too  small  compared  to  the  fire 
bed  dimensions  to  allow  attainment  of  systematic  extinction  -no 
extinction  boundaries  on  the  P  versus  q  map;  (d)  wood-crib  (i.e.. 

Class  A)  fires  also  conform  with  the  hypothesis  including  in  the 
reflashing  phenomenon,  although  the  combustion  behavior  of  crib  fires 
poses  numerous  inherent  ambiguities;  and  (e)  shredded  paper  tray  fires 
are  in  a  class  all  by  themselves  and  are  unexplained  by  the  present 
hypothesis. 


ABSTRACT 


The  mechanism  of  interaction  between  a  blast  wave  and  a  diffusion  flame 
^studied  with  the  objective  of  correlating  data  obtained  in  shocktube  tests, 
jised  on  a  scrutiny  of  the  interplay  between  physical  transport  processes  and 
cftemlcal  kinetic  reaction  processes,  a  Damkohler  number  P  is  identified  to 
represent  the  inverse  of  the  blast  strength.  Also  identified  Is  a  fire 
strength  parameter  q  which  defines  the  fuel  character.  On  a  map  of  P(q), 
regimes  in  which  extinction  is  expected  are  distinguished  from  those  in  which 
the  fire  will  sustain  the  Impact  of  the  blast.  Where  charring  fuels  are 
involved,  the  role  of  glowing  combustion  in  the  reflashing  of  an  extinguished 
flame  is  examined. 

The  hypothesis  is  tested  with  the  experimental  data  obtained  by 
Backovsky,  Martin  and  McKee  at  SRI  International  to  reach  conclusions:  (a) 
liquid  fuel  (i.e..  Class  B)  fires  are  consistent  with  the  hypothesis  (b) 
barricades  placed  upstream  of  the  fire  tend  to  produce  recirculation  of  flow 
and  as  a  result  Increase  the  flow  residence  time  to  enhance  the  stability  of 
the  fire;  (c)  the  barricade  heights  for  which  data  are  now  available  are  too 
small  compared  to  the  fire  bed  dimensions  to  allow  attainment  of  systematic 
extinction  -no  extinction  boundaries  on  the  P  versus  q  map;  (d)  wood-crib 
(I.e.,  Class  A)  fires  also  conform  with  the  hypothesis  including  in  the 
reflashing  phenomenon,  although  the  combustion  behavior  of  crib  fires  poses 
numerous  Inherent  ambiguities;  and  (e)  shredded  paper  tray  fires  are  in  a 
class  all  by  themselves  and  are  unexplained  by  the  present  hypothesis. 
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NOMENCLATURE 


A  area  of  transfer,  also  a  property  constant 

a  stoichiometric  coefficient 

B  a  property  constant 

b  stoichiometric  coefficient,  also  stick  thickness 

co  speed  of  sound 

Cp  specific  heat 

Da  Damkohler  number,  general  definition 

E  activation  energy 

f  stoichiometric  fuel/oxidant  mass  ratio 

hQ  enthalpy  of  combustion 

hp  mass  transfer  coefficient 

ko  prexponential  factor 

L  latent  heat  of  pyrolysis 

t  reaction  space  size 

m,n  orders  of  reaction 

p  pressure 

P  Damkohler  number 

q  dynamic  pressure 

q  fire  strength  parameter 

q"  fire  to  fuel  heat  flux 

R  universal  gas  constant 

T  temperature 

t*  characteristic  time 

t°  stick  burning  time 

u  velocity 


W  mass  transfer  rate 

Y  mass  fraction 

a  thermal  dlffusivity 

p  density 

y  ratio  of  gas  specific  heats 
e  nondimensional  temperature 

t  charring  solid  fire  strength  parameter,  nondimensional  preburn  time 

Subscripts 
A  fuel 

B  oxidant  (also  denoted  sometimes  by  0) 
c  char 

f  flaming 

G  glowing 

g  gas  phase 

1  initial  or  supply  state 

i  at  coordinate  *  z 

0  at  coordinate  3  zero 

0  original  solid 

p  pyrolysis 

s  surface  glowing  reaction 

•  ambient 
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INTRODUCTION 


Emanation  of  thermal  radiation  and  initiation  of  a  blast  wave  [1-4]  are 
two  of  the  early  effects  of  the  fireball  of  a  nuclear  explosion.  The  thermal 
radiation,  propagating  at  the  speed  of  light,  impinges  on  various  objects 
around  the  ground-zero.  Ignition  is  possible  within  a  critical  distance  in 
which  the  heat  fluxes  are  high  enough  and  exposures  are  long  enough.  The 
primary  fires  thus  initiated  are  established  and  made  to  grow  in  intensity  by 
a  series  of  highly  transient  physico-chemical  processes:  the  solids  continue 
to  be  transiently  heated;  combustible  gases  and  vapors  continue  to  be  formed 
by  pyrolysis  or  vaporization;  vapors  mix  with  air;  exothermic  reactions  in  the 
mixture  lead  to  evolution  of  the  flame;  and  these  flames  propagate  over  the 
burning  material . 

As  Backovsky  et  al  point  out  [3],  fire  in  this  setting  is  unlike  the 
prompt  effects  of  a  nuclear  explosion  in  that  it  continues  for  some  time  to 
take  its  own  course  of  development  and  spread  to  destroy  property  and  life. 

By  this  course  of  development,  it  is  also  perhaps  control lahle  through  the 
available  post-explosion  time,  although  its  magnitude  may  make  the  efforts  to 
control  appear  futile  in  the  wake  of  a  nuclear  attack. 

There  is  considerable  literature  on  radiant  ignition  of  idealized 
surfaces  and  some  on  temporally  distributed  exposure  and  on  nonideal  surfaces. 
There  is  extensive  research  activity  sponsored  by  the  National  Bureau  of 
Standards  [5]  on  the  topic  of  fire  growth  in  enclosures.  Fire  spread  in  the 
forest  environment  Is  a  process  relatively  well -understood  for  the  purposes  of 
control  and  damage  mitigation.  Progress  is  thus  being  made  in  gaining  an 
understanding  of  the  fire  growth  behavior  in  the  absen-.s  of  blast  wave 
effects . 

As  the  primary  fires  in  the  open  and  within  structures  continue  to  grow, 
the  blast  wave  arrives  to  perturb  them.  An  already  complicated  transient  fire 
development  process  is  now  rendered  even  more  complex  by  the  propagating 
fields  of  pressure  and  flow  brought  about  by  the  blast  wave.  The  abrupt 
adiabatic  compression  of  air  results  in  an  increase  in  temperature.  There 
also  may  be  the  tail-end  of  the  thermal  radiation  pulse.  The  growing  fire  is 
thus  subjected  to  changing  pressure  with  winds  of  varying  temperature  with  an 
externally  imposed  transient  thermal  pulse. 


Leave  alone  the  blast  damage  to  structures  for  the  moment;  enough 
knowledge  does  not  exist  on  the  manner  In  which  the  blast  winds  perturb  the 
nascent  fires.  Herein  lies  an  Important  element  of  uncertainty  in  estimating 
the  magnitude  of  the  fires  resulting  from  a  nuclear  blast.  It  Is  this  element 
which  constitutes  the  subject  of  present  research. 

The  problem  of  blast  Interaction  with  fire  was  studied  theoretically  by 
Fendell  at  TRW  [6]  and  experimentally  by  Martin,  Backovsky  and  Goodale  at  SRI, 
[1-3],  both  under  the  sponsorship  of  FEMA.  Attention  here  Is  focussed  on  the 
experimental  data  from  the  SRI  shock  tube  facility.  Data  from  both  the  shock 
tube  facility  and  the  scanty  field  tests  involving  high  explosive  simulations 
have  resulted  in  a  considerably  Improved  understanding  of  the  blast/fire 
Interaction  process.  Because  of  cost  considerations,  however,  the  available 
experiments  will  always  be  limited  In  number,  and  in  scope,  to  hamper 
deduction  of  reliable  generalizations  useful  In  FEMA's  objective  of  damage 
mitigation  and  fire  control.  Experiments  are  necessary  to  Identify  processes 
and  to  verify  theoretical  hypotheses.  The  objective  of  the  present  work  Is  to 
develop  hypotheses  of  blast/fire  Interaction  based  on,  and  to  be  verified  by, 
the  available  experiments.  Since  an  unusual  mutual  Imposition  of  physics  and 
chemistry  Is  Involved  In  the  problem  of  combustion  perturbed  by  a  blast,  an 
account  of  the  various  variables  in  a  systematic  framework  Is  desirable.  The 
research  reported  here  deals  with  such  an  account,  to  achieve  an  Improved 
consistency  among  the  available  blast/fire  interaction  data,  to  arrive  at  a 
rational  and  generallzable  interpretation  of  the  observations  to  date,  and  to 
alleviate,  at  least  in  part,  the  current  uncertainty  in  estimating  the 
magnitude  of  primary  fires  resulting  from  a  nuclear  fexplosion. 

There  are,  of  course,  other  elements  of  uncertainty  with  which  the 
decision-makers  of  emergency  management  have  to  co*pe.  These  additional  facets 
of  uncertainty  stem  from  happenings  subsequent  to  the  arrival  of  the  blast. 

The  structures  may,  and  will,  respond  to  the  blast  wave  soon.  Windows  break; 
room-filling  ensues;  loading  of  the  walls  and  floors  arises  due  to  the 
transient  overpressure  pulse;  they  collapse  by  bending  or  by  shearing;  their 
fragments  flying  around;  the  fuel  is  now  merely  piles  of  debris.  The  scene 
then  Is  probably  dominated  by  fires  in  debris,  apparently  smoldering  and 
dying,  and  yet,  apparently  eager  to  flare  up.  Winds  continue  to  blow,  and  the 
residual  thermal  radiation  Is  almost  gone.  Secondary  fires  arise  due  to 
electrical  shorts,  due  to  upsetting  of  heating  elements  In  furnaces  and  kilns 
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and  due  to  release  of  highly  flammable  fuel  liquids  and  gases  and  other 
chemicals.  Little  is  known  quantitatively  of  this  secondary  fire  scenario. 

It  is  also  in  this  stage  that  the  multiple  ^rst  effects  come  to  be  important. 
The  work  reported  here  does  not  address  the  secondary  fire  problem,  nor  the 
multiburst  problem. 

The  general  objective  of  this  work  thus  is  to  study  the  mechanisms  of 
interaction  between  a  blast  wave  and  such  a  diffusion  flame  as  expected  in 
connection  with  a  nuclear  explosion  and  as  encountered  in  simulations 
employing  shocktube  experiments  and/or  high  explosives  field  tests.  The 
specific  objective  is  to  identify  the  hypotheses  underlying  the  augmentation 
or  extinguishment  of  diffusion  flames  by  blast  waves  so  that  the  expected 
important  variables  may  be  systematically  reduced  to  nondimensional  ratios. 

The  existing  data  will  then  be  correlated  nondimensional ly  to  explain  the 
currently  available  observations  and  to  design  future  experiments. 

This  objective  is  sought  to  be  accomplished  through  the  following 
sequence  of  steps: 

(a)  systematically  identify  the  variables  Involved; 

(b)  formulate  hypothetical  mechanisms  of  flame  augmentation  or 
extinguishment; 

(c)  assess  from  these  hypotheses  the  combinations  of  variables  which 
operate  cooperatively  or  antagonistically; 

(d)  derive  from  this  assessment,  dimensionless  ratios  embodying  the 
mechanisms; 

(e)  develop  correlations  of  the  existing  data  by  means  of  these 
dimensionless  parameters; 

(f)  Identify  the  most  probable  hypotheses  of  augmentation  or 
extinguishment  which  appear  to  explain  the  existing  experimental 
observations;  and 

(g)  designate  the  specific  data  that  need  to  be  recorded  in  the  future 
blast  fire  and  shock  tube  experiments. 


APPROACH  SYNOPSIS 


Considered  here  are  fires  supported  by  hydrocarbon  liquids  and  cellulosic 
solids.  Fuels  are  characterized  by  their  thermal,  thermodynamic,  and 
combustion  properties.  The  burning  of  cellulosic  solids  takes  considerable 
time  to  become  established,  due  to  the  relatively  slow  heat  conduction. 
Furthermore  the  mechanism  of  combustible  production  is  chemical  in  nature  when 
the  fuel  is  a  cellulosic  solid  and  physical  when  it  is  a  hydrocarbon  liquid. 

As  a  result  of  these  differences,  the  burning  process  attains  a  steady  state 
quickly  If  the  fuel  is  a  liquid  and  not  so  quickly  if  it  is  a  solid.  During 
the  Initial  transient  development  phase,  thus,  the  burning  of  a  cellulosic 
fuel  array  becomes  progressively  more  intense  with  time  since  ignition.  In 
assessing  the  effects  of  any  imposed  perturbation  on  such  a  fire,  then,  the 
time  at  which  the  perturbation  is  imposed  becomes  an  important  factor.  This 
time  is  here  termed  the  'preburn  time'. 

Several  geometries  of  the  burning  surfaces  relative  to  the  flow  Induced 
by  the  blast  are  considered.  These  include:  simple  shear  flow  geometry 
( 1  -e . ,  a  flat  plate  parallel  to  the  flow);  stagnation  geometry  (i.e.,  a  plate 
normal  to  the  flow);  and  a  combined  geometry  (i.e.,  an  array  of  fuel 
cylinders  arranged  in  the  fashion  of  a  crib).  Effects  of  barriers, 
recirculation  zones  and  other  classical  flame-holding  systems  are  discussed  In 
the  light  of  these  geometries. 

The  fundamental  nature  of  burning  of  objects  is  governed  by  a  combination 
of  physics  and  chemistry  [7].  The  physical  processes  involve  preparation  of 
the  fuel  to  a  form  suitable  for  combustion  (i.e.,  vaporization,  preheating, 
e*“.),  mixing  of  the  reactants  to  form  a  combustible  mixture,  removal  of 
p<  jiucts  and  heat  from  the  flame  zone,  etc.  The  primary  chemical  process 
pertains  to  the  oxidation  reaction  between  the  fuel  and  oxygen  to  form  the 
products  and  to  release  energy.  (Production  of  combustible  gases  from  the 
fuel  (especially  if  it  were  a  solid)  may  itself  be  a  chemical  degradation 
process.)  When  the  chemical  processes  are  extremely  fast  compared  to  the 
physical  processes,  what  is  known  as  a  physically-controlled  flame  would 
result.  The  ratio  of  the  physical  characteristic  time  to  the  chemical 
characteristic  time  (or,  simply  the  chemical  rate  to  physical  rate)  is  known 
as  Damkohler  number.  Da.  When  Da  is  very  large,  a  well-established  physically 
controlled  flame  exists.  Any  changes  caused  by  the  blast  wave  In  the 
pressure,  velocity  and  the  temperature  fields  will  alter  the  relative  roles 
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played  by  chemistry  and  physics  in  the  governance  of  the  flame  behavior.  When 
the  Damkohler  number  is  made  sufficiently  low  by  the  blast  wave,  the  chemical 
processes  become  too  slow  compared  to  the  accentuated  physical  processes  and 
the  consequence  is  an  extinguishment. 

Whereas  the  Damkohler  number  plays  a  dominant  role  in  the  response  of  a 
flame  to  the  impact  of  a  blast  wave,  other  factors  may  intrude  to  modify  the 
behavior.  Barriers  and  other  recirculation  devices,  for  example,  prolong  the 
physical  residence  of  the  reactants  in  the  reaction  zone.  This  prolonged 
residence  corresponds  to  the  slowing  down  of  physical  processes.  In 
conjunction  with  a  chemical  reaction  of  a  specified  characteristic  speed,  the 
slowed  physical  processes  result  In  an  increase  in  Da.  If  this  increase  is 
sufficiently  large,  the  flame  will  become  resistant  to  blow  off.  One  ca" 
easily  write  down  the  equations  of  unsteady  conservation  of  mass,  momentum, 
energy,  and  species  for  the  boundary  layer  combustion  process  with  the 
appropriate  boundary  and  initial  conditions  [8],  An  examination  of  these 
equations,  without  actually  solving  them  readily  results  in  a  set  of 
dimensionless  criteria  useful  in  correlating  the  observed  Ignition  and 
extinction  behavior  of  the  flames. 

The  approach  used  here  Is  to  examine  the  existing  fire  extinguishment/ 
augmentation  data  in  the  light  of  different  plausible  mechanisms  of 
extinguishment,  blow-off,  rekindling,  augmentation  and  alteration  of  the 
flames  by  blast  waves.  The  excerclse  will  be  iterative  in  nature  in  that  the 
postulated  mechanisms  attempt  to  explain  and  correlate  the  data  while  the  data 
and  observations  direct  towards  identifying  the  plausible  mechanisms. 


PREPARATORIES 


What  is  a  Blast  Wave? 

The  generation  of  a  pressure  wave  due  to  rapid  energy  release  is  called 
an  explosion  [4].  When  a  0.5  kg  sphere  of  TNT  (whose  volume  is  about  3xl0-4 
m^)  is  activated,  about  3x10"!  m3  of  combustion  gases  at  about  4  000  K  and  100 
000  atm  are  rapidly  produced.  These  high  pressure,  high  temperature  gases 
isentropically  expand  to  'shock  up'  the  surrounding  compressible  gas,  namely 
the  air.  A  spatial,  (generally  moving),  discontinuity  known  as  a  shock  front 
is  thereby  formed,  across  which  pressure,  temperature  and  density  exhibit 
abrupt  jumps.  Passage  of  the  shock  front  through  an  otherwise  still  air 
causes  air  to  move,  thus  producing  a  particle  velocity.  The  speed  of  winds 
thus  induced  is  frequently  expressed  as  a  dynamic  pressure.  The  quantity  of, 
and  the  rapidity  with  which,  energy  is  released  by  the  source  determines  the 
intensity  of  the  shock  front,  i.e.,  shock  velocity,  peak  over  pressure, 
dynamic  pressure,  temperature  rise,  etc.  As  the  shock  wave  encounters 
obstacles  in  the  course  of  its  propagation,  it  is  partially  reflected  (to 
augment  the  overpressure)  and  partially  diffracted  around  the  obstacle 
(imparting  some  of  its  energy  to  the  obstacle).  One  of  the  consequences  of 
this  diffraction  is  the  evolution  of  a  recirculating  (or  vortex)  flow  on  the 
downstream  side  of  the  obstacle. 

Books  on  gas  dynamics  present  a  quantitative  description  of  the  shock 
propagation  in  various  media.  Glasstone  summarizes  the  relevant  relations  for 
nuclear  explosions  in  Reference  4. 

Near  time  equal  to  zero,  the  fire  ball  produced  by  an  explosion  is  such 
that  pressure  is  high  near  its  edge  and  is  progressively  smaller  as  its  origin 
is  approached.  Upon  expansion  of  the  fire  ball,  the  overpressure  starts  to 
decay  so  that  after  a  certain  time,  the  pressure  behind  the  shock  falls  to  the 
ambient  pressure.  At  subsequent  times,  a  partial  vacuum  is  formed  behind  the 
shock  front.  This  is  called  the  negative  phase.  While  winds  blow  in  the 
direction  of  the  shock  front  propagation  during  the  duration  of  the  positive 
phase  ,  they  blow  towards  the  source  during  the  negative  phase.  The  steady 
decay  of  the  peak  overpressure  with  time  (or  distance  traversed  by  the  shock) 
is  presented  in  Fig.  3.72  of  Ref.  4.  As  the  peak  overpressure,  the  positive 
phase  duration  as  well  as  the  negative  phase  intensity  and  duration  are 
dependent  upon  the  distance  travelled  by  the  shockwave.  While  peak 


overpressure  is  an  indicator  of  blast  intensity,  dynamic  pressure  (pu^/2)  is 
an  indicator  of  the  wind  velocity  u.  (p  is  air  density  behind  the  shock). 
Peak  dynamic  pressure  also  decays  with  increasing  distance  from  the  source. 

The  peak  overpressure  decay  (when  coupled  with  appropriate  scaling  laws 
for  the  source  strength  and  height  of  burst*  and  with  the  following  gas 
dynamic  relations  from  Reference  4)  leads  to  estimation  of  a  variety  of  blast 
characteristics  of  relevance  in  the  current  study. 


^reflected  2p  + 

where  Po,  po  and  co  are  respectively  the  ambient  air  pressure,  density  and 
speed  of  sound,  p  is  peak  overpressure,  y  is  ratio  of  specific  heats  («1.4 
for  air),  u  is  velocity,  p  is  density,  q  is  dynamic  pressure,  and  Preflected 
is  reflected  overpressure  (for  normal  incidence  on  a  flat  surface). 

The  energy  released  at  the  explosion  source  is  distributed  among  a 
variety  of  modes.  These  modes  include,  first,  the  hot  products  of  explosion 
thermally  radiate  to  the  surroundings  in  a  characteristic  pulse  [4],  Second, 
the  expanding  hot  gas  performs  mechanical  work  on  the  air  to  produce  and 
propagate  the  shock  wave  (which  itself  imparts  part  of  its  energy  to  raise  the 
kinetic  energy  of  air  by  producing  winds  which  would  eventually  be 

^  The  primary  scaling  consideration  is  that,  for  a  given  peak  overpressure  or 
dynamic  pressure,  the  height  of  burst,  distance  from  ground  zero,  time  of 
blast  arrival,  and  the  duration  of  positive  phase,  are  all  proportional  to  the 
cube-root  of  the  weapon-yield  (which  is  generally  expressed  in  units  of 
kiloton  TNT  equivalence). 


dissipated).  Furthermore,  the  isentropic  compression  of  air  results  in  a 
temperature  rise.  The  sensible  energy  increase  associated  with  the 
temperature  rise  is  irrecoverably  lost.  A  part  of  the  energy  involved  in  the 
winds  and  static  pressure  will  be  imparted  to  fragments  of  debris  formed  from 
the  break  up  of  objects  in  the  path  of  the  blast  wave;  these  fragments  gain 
kinetic  and  potential  energies  while  being  physically  displaced  through 
significant  distances.  The  bursting  of  objects,  loaded  by  the  blast  to  form 
the  fracpnents  is  itself  a  process  in  which  energy  is  absorbed  to  result  in 
strain  and  rupture. 

What  is  a  Fire? 

Strongly  exothermic  oxidations  constitute  the  bulk  of  combustion 
reactions.  These  reactions  may  occur  in  the  gas  phase,  at  the  solid  surface 
or  within  a  porous  solid  to  result  respectively  in  flaming,  glowing  and 
smoldering.  Heating  of  the  fuel,  first  by  an  external  means  and  then  by  the 
energy  release  by  combustion,  to  produce  combustible  vapors,  is  a  predominant 
feature  of  the  combustion  of  solids  and  liquids.  While  liquid  fuel  fires  are 
known  as  Class  B  fires,  those  of  wood-like  charring  solids  are  known  as  Class 
A  fires.  A  fire,  in  a  broad  sense,  is  a  collection  of  flames  generally 
accompanied,  in  Class  A  systems,  by  glowing  as  well  as  smoldering.  If  the 
fuel  were  supplied  independently  in  gaseous  state,  the  flame  character  will 
be,  by  and  large,  externally  controlled.  If  the  fuel  otherwise  is  supplied  as 
a  liquid  or  solid,  the  flame  feeds  some  of  its  energy  back  to  the  fuel  source 
to  produce  combustibles  at  such  a  rate  as  to  maintain  a  self-imposed 
pseudo-equilibrium  in  its  characteristics.  The  workings  of  the  flame  in  both 
these  cases  involve  such  component  processes  as  mixing  of  the  fuel  vapor  with 
air,  preheating  the  mixture,  supporting  the  combustion  chemical  reaction,  and 
dispersing  the  energy  and  product  species  released  in  the  reaction.  In  a 
flame,  these  processes  are  generally  consecutive  rather  than  simultaneous.  As 
a  consequence,  the  slowest  of  these  process  steps  governs  the  overall  behavior 
of  the  combustion. 

Generally,  combustion  in  fires  is  governed  by  the  physics  of  fuel 
generation,  mixing,  and  heat  transfer  rather  than  by  the  relatively  faster 
chemical  kinetics.  In  flames  near  extinction,  however,  the  chemical  reaction 
kinetics  are  slowed  down  to  a  speed  closer  to  that  of  the  physical  transport 
pheomena.  This  issue  will  be  further  discussed  in  a  later  section. 


Ignition  of  liquid  fuel  fires  involves  preheating  the  fuel  to  produce  a 
sufficient  quantity  of  vapor  which,  when  mixed  with  air  under  the  specific 
fluid  dynamic  mixing  conditions,  yields  a  mixture  above  the  lean  ignitibility 
limit.  The  introduction  of  an  ignition  source  (a  pilot  flame,  spark,  heated 
wire  or  an  incandescent  particle)  then  produces  a  flame  which  propagates  over 
the  entire  available  liquid  surface.  This  flame  will  soon  establish  itself 
over  the  liquid  fuel,  first  going  through  a  short  transient  period  in  which 
the  liquid  is  heated  to  steady  state,  and  then  in  a  steady  state.  The  liquid 
surface  is  usually  not  hotter  than  the  saturation  temperature  corresponding  to 
the  ambient  pressure.  For  most  of  the  common  liquid  fuels,  this  boiling-point 
temperature  is  relatively  low.  If  hot  metal  rims,  incandescent  particles  or 
other  ignition  sources  are  absent,  then  a  blown-out  flame  remains  blown-out 
even  if  the  mixture  conditions  are  suitable.  Note,  however,  that  a  flame  may 
either  be  blown  away  by  a  blast  impingement,  or  be  simply  displaced 
down-stream.  A  displaced  flame  can  reflash  to  the  fuel  source  as  vapors 
continue  to  be  produced  and  as  the  blast  effects  subside. 

Consider,  as  shown  in  Fig.  1,  a  steadily  burning  pool  of  liquid  fuel. 

Fig.  la  indicates  the  pool  burning  in  quiescent  air  while  lb  shows  one  burning 
in  flowing  air.  Radiative  and  convective  heat  transfer  from  the  flame  to  the 
fuel  surface  causes  evaporation  to  provide  combustible  vapors  to  the  flame. 

Air  (oxygen)  is  transported  to  the  flame  from  the  surroundings  to  react  with 
the  fuel.  All  the  reaction  products  and  most  of  the  heat  produced  are 
convected  out  of  the  flame;  part  of  the  energy  is  also  radiated  away. 

Typical  composition  and  temperature  profiles  are  indicated  in  Fig.  2. 

The  gas  phase  combustion  chemistry  is  assumed  to  be  infinitely  fast  in  Fig.  2a 
compared  to  the  physical  processes  of  mixing  and  transfer  of  energy,  species 
and  momentum.  Under  this  assumption  the  flame  is  a  plane  of  zero  thickness  at 
which  the  reactants  undergo  combustion  as  fast  as  they  are  transported  to  it. 
The  dominant  feature  of  this  sort  of  infinite  rate  chemistry  is  that  the 
oxygen  and  fuel  are  prevented  from  penetrating  through  the  flame.  This  is 
known  as  an  ideal  diffusion  flame.  Its  characteristics  are  obtainable  through 
an  analysis  of  the  associated  fluid  mechanics  problem  without  regard  to  the 
chemical  reaction  kinetic  detail.  A  further  elaboration  of  this  will  be  given 
in  a  later  section. 
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If  one  tampers  with  the  problem  of  Fig.  2a  either  to  reduce  the  chemical 
reaction  rate  or  to  enhance  the  physical  transport  rates,  the  flame  sheet 
ceases  to  be  infinitesimally  thin,  reactants  may  now  coexist  and  the  flame 
temperature  is  lowered  as  depicted  in  Fig.  2b.  Upon  thus  sufficiently 
broadening  due  to  the  vigorous  transfer  mechanisms,  the  flame  would  altogether 
be  blown  out.  This  extinguishment  process  will  be  further  discussed  later 
since  it  has  a  bearing  on  blast/fire  interaction. 

Class  A  Fires: 

Wood  is  a  member  of  cellulosic  solid  fuel  family  which  undergoes  charring 
combustion.  Production  of  combustible  gases  (and  vapors)  upon  heating  wood  is 
a  chemical  degradation  process.  Also  known  as  pyrolysis  or  charring,  this 
mode  of  thermal  response  of  the  solid  fuel  is  in  contrast  with  the  physical 
(thermodynamic)  evaporation  of  a  liquid  fuel. 

As  the  combustible  gases  produced  in  the  pyrolysis  exit  from  the  solid  to 
the  gas  phase,  they  form  a  mixture  with  air  in  the  boundary  layer  adjacent  to 
the  solid.  This  mixture  may  be  ignited  either  spontaneously  or  by  a  pilot 
source.  Even  after  such  an  ignition,  continued  burning  may  not  be  possible 
without  an  external  heat  flux  to  the  solid.  The  reason  for  this  is  that 
energy  is  continuously  conducted  away  into  the  interior  of  the  solid  and 
radiated  out  to  the  surroundings;  a  feeble  nascent  flame  can  not  quite  provide 
sufficient  energy  for  these  conductive  and  radiative  drains  as  well  as  the 
energy  required  for  the  continuation  of  pyrolysis.  Upon  continued  external 
heating,  however,  a  stage  soon  arrives  at  which  self-sustained  burning  is 
possible. 

When  the  cellulosic  fuel  described  above  is  a  pile  of  sticks  of  wood, 
mutual  interactions  among  the  different  fuel  elements  (i.e.,  sticks)  in 
conserving  energy  would  tend  to  aid  sustained  burning.  The  longer  the 
burning,  the  more  established  the  fire  would  become. 

It  is  Important  to  note  that  the  surface  temperature  of  burning  wood  is 
significantly  higher  than  that  of  a  typical  burning  liquid  fuel.  If  the  flame 
were  to  be  blown  away  momentarily  for  whatever  reason,  the  thermal  inertia  of 
the  solid  sustains  the  pyrolysis  and  the  hot  surface  is  capable  of  piloting  a 
relgnltlon.  Additionally,  if  the  flame  were  to  be  blown  out  or  if  it  were  so 
defective  in  structure  as  to  permit  penetration  of  oxygen  to  the  hot  surface, 
glowing  combustion  may  ensue. 


The  role  of  glowing  combustion  will  be  further  discussed  elsewhere  in 
this  report.  Suffice  it  to  note  that  glowing  is  expected  to  present  a 
substantial  energy  source  at  the  surface.  Enough  of  the  energy  produced  in 
this  surface-glowing  will  be  transmitted  to  the  solid  to  keep  the  pyrolysis 
continued  and  to  thus  keep  combustible  supply  to  the  gas  phase  continued.  The 
glowing  char  surface  is  generally  capable  of  reigniting  the  gas  mixture  thus 
formed  following  the  blowing  away  of  the  original  flame.  Herein  lies  a  major 
difference  between  the  burning  of  solids  and  liquids. 

Control  by  Physics  or  Chemistry? 

As  pointed  out  above,  the  burning  of  any  fuel  involves  both  physics  and 
chemistry.  The  physical  processes  involved  are  heat,  mass  and  momentum 
transport,  thermal  radiation  and  thermodynamics.  The  oxidative  process  is 
chemical.  To  delineate  the  relative  roles  of  these  physical  and  chemical 
processes,  consider  two  infinite  parallel  plates  separated  by  a  distance  l  as 
shown  in  Fig.  3.  Let  one  of  these  plates  be  the  source  of  fuel  vapor  A  while 
the  other  is  the  source  of  the  oxidant  B.  Let  YAq  and  YBjt  be  the  mass 
fractions  of  A  and  B  near  the  sources.  As  A  and  B  diffuse  into  the  space 
between  the  plates,  chemical  reaction  will  ensue. 

The  physical  rate  of  transfer  of  the  reactant  species  from  the  source  to 
the  reaction  space  occurs  in  proportion  to  the  difference  in  mass  fraction, 
the  proportionality  constant  being  the  mass  transfer  coefficient  which  is  a 
function  of  the  flow  characteristics  of  the  reacting  medium,  the  diffusion 
coefficient  and  density  of  the  species  through  the  reacting  mixture  and  the 
characteristic  dimension  of  the  transfer  surface.  Focussing  attention  on  the 
fuel  species  for  convenience, 

«A  -  HDaA(YA0-Ya)  (1) 

where  WA  is  the  mass  of  species  A  transferred  per  unit  time  (kg/s),  h[)A  is  the 
mass  transfer  coefficient  (kg/m^s),  A  is  the  area  of  plate  (m^),  YAq  is  the 
species  A  mass  fraction  at  the  source  and  Ya  is  the  same  in  the  reaction 
space.  The  mass  transfer  rate  for  the  oxidant  can  similarly  be  written  as 

Wb  »  hDBA(YBjrYB)  •  (2) 

By  stoichiometry  aA  +  bB  -*■  products,  the  rates  are  related  by 
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d6  ■  <$>  SA  (3) 

The  chemical  kinetic  rate  is  generally  of  the  form 

WA  =  kQA  At  YAnYBm  exp(-E/RT)  (4) 

where  kg,  n,  m  and  E  are  the  reaction  kinetic  parameters  known  respectively  as 
the  prexponential  (collision)  factor,  order  of  the  reaction  with  respect  to 
fuel  A  and  to  oxidant  B  and  activation  energy.  R  is  the  universal  gas 

constant.  T  is  the  reaction  space  temperature  assumed  here  to  be  a  prescribed 
constant.  The  above  four  equations  can  now  be  solved  for  the  four  unknowns: 
the  reaction  space  concentrations  YA  and  Yg  and  the  rates  of  consumption  foA 

and  fog.  Specifically, 

VA  "  rA0  '  (HH^)  (6) 


’»  *  v  •  <*>  ‘fe  ^ 


‘  [knaAJ  exp  (j*)^ 


For  given  values  of  n,  m,  YA0,  YB*,  (hoA/koA  *  exp(-E/RT),  (b/a)  and  (hDA/hoB). 
Eq.  (7)  can  be  solved  for  and  then  Eqs.  (5)  and  (6)  give  Yft  and  Yg. 

The  general  solution  is  tedious  and  obscuring.  In  order  to  elucidate  the 
general  features  of  the  solution,  consider  n  =  1  and  m  =  0.  Then 


YA  =  yA0  TRTa 


SA  *  ^AO  TO  hDA*  ■  VA0  TO  k0AM 


where 


kn.Ai.  exp(-E/RT) 


Consider  now  two  extreme  cases  when  Da  >>  1  and  Da  «  1.  Da  is  much  greater 
than  unity  when  the  chemical  reactions  occur  very  rapidly  compared  to 
physical  transfer  of  species.  Equations  (8)  and  (9)  then  reduce  to  Ya  - 
(Y/^q/D)  *  0(s1nce  the  mass  fraction  can  at  best  be  equal  to  unity)  and 

&A  "  *A0*  Since  the  chemical  kinetics  are  infinitely  fast,  the 

comparatively  slower  transport  physics  governs  the  global  characteristics 
of  the  problem.  That  is,  the  chemical  reaction  consumes  the  reactants 
wherever  the  stoichiometric  requirements  are  fulfilled  as  soon  as  they  are 
brought  into  the  reaction  zone.  Since  the  transport  rates  are  low,  the  driving 
forces  for  transfer  tend  to  be  maximized  as  reflected  by  the  strong  spatial 
gradients  (of  the  species  mass  fraction  in  this  illustration).  The  so-called 
'sudden  death  kinetics'  culminate  in  the  reaction  at  a  distinct  site  to  which 
the  fuel  and  oxidant  are  supplied  in  stoichiometric  proportions. 

Said  to  be  physically  (or  diffusionally)  controlled  combustion,  this 
limiting  behavior  is  depicted  in  Fig.  3.  The  flame  located  at  station  marked  1 
is  an  infinitesimally  thin  imaginary  sheet.  This  sort  of  an  ideal  diffusion 
flame  arises  whenever  high  temperatures  and  pressures,  large  reactor 
dimensions,  low  fluid  velocities  and  low  species  diffusi vities  are  encountered. 
If  Da  were  not  infinity  but  still  large,  the  fuel  and  oxidant  may  not  coexist 
at  any  station;  as  a  result,  the  flame  ceases  to  be  a  sheet  of  zero  thickness. 
Said  differently,  finite  rate  chemistry  broadens  the  flame  (between  stations 
labelled  2  and  3  in  Fig.  3).  Chemistry  in  these  nonideal  diffusion  flames  can 
not  be  altogether  ignored  in  attempts  to  arrive  at  the  combustion 
characteristics. 

Consider  now  the  second  limiting  case  when  Da  «  1.  This  situation 
signifies  that  the  chemical  reaction  rate  is  extremely  small  compared  to  the 
physical  transport  rates.  Low  temperatures  and  pressure,  small  physical 
dimensions  of  the  reactor,  greatly  enhanced  flow  and  mixing  conditions  are 
conducive  to  this  situation.  Setting  D  =  0,  Eqs.  (8)  and  (9)  yield  Ya  -  Y/vo 
and  6^  «  k^At  ^A0YBt  exP( -E/RT) .  The  relatively  slower  chemical  kinetics 
govern  the  global  characteristics  of  this  limiting  combustion  while  the 
relatively  fast  fluid  dynamic,  thermodynamic  and  transfer  processes  become 
Irrelevant  simply  by  washing  out  the  spatial  gradients.  A  distinct  reaction 
zone  (l.e.,  flame)  does  not  exist.  This  kinetically  controlled  oxidation 
process  is  also  depicted  in  Fig.  3  through  the  spatially  uniform  composition 
profiles. 


Comments 


(a)  In  the  preceding  discussion,  we  have  deliberately  kept  the  reactor 
temperature  fixed  to  develop  the  concept  of  chemical  and  physical  control  of  the 
combustion  process.  A  relaxation  of  this  artificial  restraint  leads  to  an 
interesting  explanation  of  the  ignition  and  extinction  phenomena  relevant  to  the 
blast/fire  interaction.  We  will  carry  this  relaxation  in  the  following 
chapter. 

(b)  The  general  descriptive  results  of  the  above  analysis  are  valid  for 
all  values  of  n  and  m,  the  orders  of  the  reaction  with  respect  to  the  fuel  and 
the  oxidant. 

(c)  The  nondimensional  parameter  Da  is  known  as  Damkohler  number.  In  the 
preceding  discussion,  it  is  defined  as  the  ratio  of  characteristic  chemical 
reaction  rate  to  the  physical  supply  rate  of  species  A.  It  could  also  have  been 
defined  on  the  basis  of  species  B,  product  species,  or  energy.  These  various 
definitions  will  be  related  through  the  stoichiometry.  The  physical  rate 
appearing  in  the  denominator  of  any  of  these  definitions  may  be  based  on 
molecular  diffusion,  convective  transport,  turbulent  diffusion,  or,  simply,  the 
flow  rate.  The  Da  values  obtained  by  these  various  physical  rates  are  mutually 
related  through  the  thermo-fluid  mechanical  relations  among  the  different 
physical  transport  quantities.  No  matter  the  specific  details  of  definition, 
the  Damkohler  number  is  always  a  ratio  of  the  chemical  speed  to  physical  speed, 
always  indicating  the  relative  magnitudes  of  chemical  to  physical  rates. 


EXTINCTION  AND  IGNITION 


Of  concern  in  this  work  is  the  behavior  of  fires  over  liquid  fuel  pools 
and  wood  fuel  cribs  upon  which  a  blast  wave  is  imposed.  We  will  discuss  in 
this  chapter  the  interaction  of  a  blast  with  liquid  fuel  pool  fire  and  then 
extend  the  treatment  to  a  wood  crib  fire  which  involves  pyrolysis  and  glowing 
in  addition  to  gas  phase  flames. 

Pool  Fire: 

Consider  a  steadily  burning  pool  fire  shown  in  Fig.  la.  As  a  blast  wave 
is  impinged  on  this  fire  with  its  attendant  perturbations  in  wind,  pressure  and 
temperature,  the  flame  will  be  disturbed.  Such  a  disturbance  may  manifest  in 
one  or  more  of  the  following  physico-chemical  phenomena. 

(a)  Annihilation  of  spacial  gradients  of  species,  temperature  and  velocity 
by  the  increased  molecular  and  turbulent  transport  can  lead  to  intense  mixing 
and  dilution  of  the  reactants.  Also  brought  about  is  thermal  dilution  and  a 
consequent  reduction  in  the  temperature  of  the  reaction  space.  The  net  overall 
effect  is  to  reduce  the  Damkohler  number. 

(b)  Energy  feedback  from  the  flame  to  the  condensed  phase  fuel  will  be 
reduced  due  to  physical  displacement  or  deformation  of  the  flame.  This 
reduction  would  result  in  a  decrease  in  fuel  vapor  supply  to  the  gas  phase.  The 
consequent  decrease  in  fuel  vapor  concentration  may  reduce  the  chemical  kinetic 
rate  to  yield  a  diminished  Damkohler  number.  A  displaced  or  deformed  flame 
could  also  be  cause  for  a  reduction  of  the  temperature  of  gas  phase  near  the 
fuel  bed,  thus  accentuating  the  decrease  in  Da. 

(c)  Energy  feedback  to  the  fuel  bed  could  be  augmented  by  the  wind  bringing 
the  flame  closer  to  the  surface.  A  higher  wind  velocity  would  then  reduce  the 
boundary  layer  thickness  and  flame  stand-off  distance  to  increase  the  heat 
transfer.  This  effect  can  be  expected  only  when  the  wind  is  feeble.  The  overall 
result  is  an  augmentation  of  the  fire  intensity  by  the  imposition  of  the  blast. 
This  effect  is  contrary  to  that  discussed  in  item  (b)  above. 

(d)  Energy  feedback  to  the  fuel  bed  will  be  enhanced  due  to  the  flame 
holding  in  the  reci rculatory  zones  of  a  baffled-barrier-flow.  The  barrier  could 
be  the  rim  of  a  fuel  pan  or  simply  a  step  located  upstream  of  the  pool . 

(e)  The  fuel  bed  may  be  mechanically  disturbed  and  dispersed  by  the  blast 
wave.  Such  dismemberment,  dispersion  and  redistribution  of  the  fuel  bed  may. 
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under  some  conditions,  increase  the  fire  intensity  by  transforming  the  bulk  fuel 
confined  in  a  tray  into  an  unconfined  spray.  Under  other  conditions,  however, 
fragmentation  of  the  fuel  bed  might  aid  to  dissipate  the  energy  content  of  the 
fuel  in  the  tray  to  an  ineffectually  low  average  level.  In  the  first  situation, 
the  dispersion  to  form  a  spray  would  enhance  the  transfer  surface  area  available 
for  fuel  vaporization.  Energy  interaction  among  fragments  of  fuel  is 
significant.  In  the  second  situation,  the  same  increase  in  transfer  area  will 
make  the  fuel  fragments  lose  energy  to  the  air  flow  because  the  mutual 
interactions  are  minimal.  Although  significantly  strong  blasts  are  expected  to 
produce  the  second  situation,  precise  delineation  of  the  rare  conditions 
required  to  accomplish  this  is  not  easy.  It  appears  that  one  can  expect  the 
first  situation  almost  always  in  practical  blast/fire  interaction  conditions. 

(f)  Pressure  change  will  result  in  a  shift  in  combustion  chemical 
kinetics.  By  and  large,  combustion  of  hydrocarbons  in  air  is  approximately  a 
second  order  reaction.  As  a  result,  the  kinetic  rate  is  proportional  to  the 
square  of  pressure,  so  that  the  Damkohler  number  increases  with  increasing 
pressure.  The  increase  in  pressure  also  alters  the  fluid  dynamics  to  increase 
the  coefficients  of  heat  and  mass  transfer,  so  that  the  Damkohler  number  is 
reduced.  The  net  effect  of  these  two  opposing  effects  of  increased  pressure  can 
not  be  predicted  without  a  detailed  study  of  the  underlying  basic  coupled 
phenomena.  Note,  however,  that  the  pressure  rise  associated  with  a  blast  wave 
is  temporally  variant.  As  a  result,  arguments  based  on  static  imposition  of  a 
pressure  rise  might  become  invalid  in  the  dynamic  behavior  of  the  blast-impacted 
flame. 

(g)  The  shockwave  also  brings  with  it  a  temperature  rise  due  to  isentropic 
compression  of  air.  This,  too,  is  a  transient  phenomenon.  It  can  be  expected 
to  have  some  effect  on  the  flame  especially  in  its  chemical  kinetic  aspects. 

(h)  In  all  practical  situations  of  blast  wave  generation  by  the  explosion 
of  a  weapon,  a  thermal  radiation  pulse  is  involved  [4],  The  effect  of  the 
residual  radiant  pulse  is  to  contribute  to  vaporization  of  the  fuel  bed  even 
when  the  energy  feedback  is  mitigated  from  the  displaced  or  disappeared  flames. 
Even  more  important  is  the  thermal  radiation  pulse  from  subsequent  weapon  bursts 
in  a  multiburst  scenario.  The  effect  of  a  radiation  pulse  is  altogether  ignored 
in  the  present  work. 

By  and  large,  most  of  the  effects  discussed  above  can  be  condensed  to  the 
following  three  global  alterations  brought  about  by  the  impingement  of  a  blast 
on  a  flame:  first,  the  gas  phase  reaction  domain  (which  was  occupied  by  the 
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flame  prior  to  blast  arrival)  suffers  excessive  thermal  dilution;  second,  the 
fuel  vapor  concentration  in  this  domain  is  diminished;  and  third,  oxygen 
concentration  of  the  domain  is  increased.  The  first  and  second  items  dominate 
in  determining  the  fate  of  the  flaming.  The  flame  is  generally  blown  out 
completely  or  displaced  to  a  downstream  location.  The  hot  fuel  layer  in  the 
condensed  phase  loses  heat  to  the  gas  phase  by  convection  and  to  the  interior 
of  the  fuel  bed  by  conduction/convection.  The  temperature  of  both  the  gas 
phase  and  the  fuel  bed  decrease  gradually  and  monotonically  with  time;  so  also 
does  the  gas  phase  fuel  concentration.  If  the  blast  effects  subside  before  the 
thermal  and  species  dilution  reach  critical  limits,  the  fire  may  reflash  or 
reignite.  To  find  these  critical  limits  is  the  objective  of  the  following 
simple  analysis. 


Analysis: 

The  interaction  process  can  be  studied  rigorously  by  solving  the  equations 
of  conservation  of  total  mass,  momentum,  energy  and  species  under  the 
constraints  of  known  boundary  and  initial  conditions  along  with  the 
supplementary  relations  pertaining  to  turbulence,  saturation  thermodynamics, 
ideal  gas  law,  laws  of  gas  mixtures,  transport  laws,  chemical  reaction 
kinetics,  energetics  and  equilibria  and  stoichiometry.  These  general 
conservation  equations  are  available  in  combustion  textbooks.  All  combustion 
problems  are  but  special  cases  described  by  these  equations.  With  respect  to 
ignition  and  extinction,  one  seeks  solutions  for  the  problem  in  which  energy 
and  species  concentrations  are  so  disturbed  by  the  flow,  to  cause  large  and 
abrupt  temperature  rise  (ignition)  or  decrease  (extinction).  Due  to  the 
inherent  nonlinearities  in  the  governing  equations  only  numerical  solutions  are 
generally  possible  although  some  asymptotic  solutions  of  the  limiting  cases  are 
in  existence.  In  the  following,  a  simple  algebraic  analysis  is  presented  to 
obtain  a  closed-form  solution  of  the  reacting  flow  problem.  This  solution 
clearly  depicts  the  ignition  and  extinction  phenomena. 

Suppose  that  Y/^  is  the  mass  fraction  of  the  fuel  at  the  supply  boundary 
into  the  reacting  gas  phase  in  which  the  mass  fraction  is  Ya.  If  the  average 
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reaction  rate  is  (-W^  )  kg/m  s,  the  characteristic  reaction  time,  t*.  is  of 

the  order 

p^Ai"^A^ 

t* - -  (11) 
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The  reaction  rate  may  be  written  as 


-WA’ 1 '  -  k  P  Ya 

(12) 

k  h  kg  exp(-E/ET) 

(13) 

where  the  reaction  is  taken  to  be  of  first  order  for  illustrative  purpose,  k  is 

the  specific  reaction  rate  constant  {units  1/sec),  kg  is  the  preexponential 

collision  factor  (units  1/sec),  E  is  activation  energy  (J/kgmole),  R  is  the 
universal  gas  constant  (J/kgmole  K);  and  T  is  the  characteristic  temperature  (K) 

of  the  reaction  system,  p  is  the  density  reacting  mixture,  (kg/m^),  a  function 

of  the  pressure,  temperature  and  composition  of  the  reacting  gas  phase. 

Combining  Eqs.  (11)  and  (12)  the  fractional  destruction  of  fuel  species  by 

combustion  is  given  by 


YAi"YA  kt* 
YAi  '1+kt* 


(14) 


t*  is  to  be  equal  to  the  time  provided  by  physics.  As  such,  it  may  be  expressed 

in  various  ways.  It  can  be  the  reaction  volume  (V)  divided  by  the  volumetric  flow 

rate  of  reactants  (tf).  It  can  also  be  the  ratio  of  reaction  path  length  (t)  to 
the  flow  velocity  (u).  Equally  well,  it  can  be  the  length  i  squared  over  the 

diffusion  coefficient.  Other  ways  of  defining  the  time  are  also  suitable.  Thus, 

Eq.  (14)  is  merely  the  species  conservation  principle  in  which  the  term  kt* 

corresponds  to  the  ratio  of  chemical  reaction  rate  to  a  physical  rate  (i.e., 

diffusion  rate,  convection  rate,  etc.)  The  energy  conservation,  similarly,  is 

given  by 


VAi-',A 


<WT-V 


(15) 


Where  Cpg  is  the  gas  phase  specific  thermal  capacity  (J/kg  K),  T,  is  the  source, 
reference,  or  initial  temperature  and  hc  is  enthalpy  of  combustion  (J/kg  fuei). 
The  origin  of  Eq.  (15)  is,  in  fact,  a  combination  of  the  species  A  and  energy 
conservation  equations  between  which  the  nonlinear  reaction  terms  can  be 
eliminated  by  Schvab-Zeldovitch  transformation  [7,8]  involving  the  definition 
of  a  composite  variable  [hcYA  +  Cpg  T],  conserved  under  the  usual  convective 
diffusion  constraints. 

From  Eqs.  (14)  and  (15)  one  can  now  obtain 

<T-V  kt* 

*A1  hc  ‘  l+M* 
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so  that  with  Eq.  (13), 


P  = 


e1/0(0-ei) 

q-(9-ei) 


(16) 


where,  with  t  =  fc/u  =  the  time  taken  by  the  flow  of  velocity  u  to  traverse  a 
lengh  I , 

e  =  RT/E  ei  =  RT-j /E 

P  i  k0t/u  q  H  RhcVA./ECpg 

The  quantity  P  is  a  Damkohler  number,  q  is  indicative  of  the  energetic 
strength  of  the  reaction.  0-j  is  the  initial  or  boundary  temperature  normalized 
with  the  reaction  activation  temperature.  Whereas  q  gives  a  measure  of  energy 
release  quantity,  the  nondimensional  activation  energy  ( 1/e ^ )  gives  a  measure 
of  the  temperature-dependency  of  the  combustion  reaction. 

In  Figure  4,  a  sketch  of  Eq.  (16)  is  given  to  demonstrate  the  general 
character  of  0  dependence  of  P  for  any  given  values  of  q  and  ©i.  The  curve  is 
generally  s-shaped.  The  lowest  branch  of  this  curve,  insensitive  to  both  P  and 
q,  corresponds  to  the  extremely  slow  oxidation  reactions  and  is  known  as  the 
'frozen'  branch.  As  one  increases  the  chemical  reaction  rate,  P  increases  and 
minute  temperature  rise  will  result.  Upon  thus  reaching  the  state  identified  by 
I,  Eq.  (16)  says  that  the  reaction  system  has  to  follow  the  course  represented 
by  the  broken  line  in  Fig.  4.  This  course,  however,  is  physically  forbidden 
since  it  implies  that  a  decrease  in  reaction  rate  or  an  increase  in  physical 
dispersion  rate  (both  of  which  make  P  fall)  would  increase  the  reaction 
temperature.  The  only  physically  allowed  recourse  for  the  reaction  process  then 
is  to  experience  an  abrupt  jump  from  state  I  to  state  I*  on  the  high  temperature 
branch  of  the  solution  which  stands  for  the  equilibrium  combustion  resulting  in 
a  flame.  This  abrupt  jump  is  known  as  ignition. 

Consider  now  a  flame  already  existent  somewhere  on  the  equilibrium  branch 
of  the  solution.  As  this  flame  is  so  perturbed  by  external  action  as  to  reduce 
P  either  by  decreasing  the  numerator  involving  the  characteristic  chemical 
kinetic  rate  or  by  increasing  the  denominator  involving  the  characteristic 
physical  transport  rate,  the  flame  temperature  gradually  diminishes  until  state 
E  is  reached.  From  state  E,  the  system  is  asked  to  enter  the  forbidden  branch. 
As  a  result,  it  jumps  abruptly  from  E  to  E'  on  the  frozen  branch,  this  jump 
being  associated  with  what  is  known  as  extinction. 


A  branch  of  applied  mathematics  known  as  catastrophe  theory  evolved  over 
the  last  two  decades  to  describe  the  abrupt  jump  processes  of  this  sort  as  they 
arise  in  the  fields  of  chemical  engineering,  biology,  elasticity,  and  disaster 
management.  Problems  involving  jump  processes  in  these  fields  are  generally 
not  amanable  to  solutions  by  the  classical  methods  of  mathematics  of  continuous 
functions. 

Having  thus  described  the  general  nature  of  the  solution,  we  present  in 
Figs.  5-8  the  parametric  dependency  of  the  solution.  Figure  5  indicates  that 
the  hysterisis  between  II 'EE'  discussed  above  is  not  apparent  when  the 
activation  energy  is  low,  for  then  the  temperature-dependence  of  the  reaction 
rate  is  rather  weak.  (Recall  that  6,  e-j,  and  q  all  become  large  as  E  becomes 
small).  The  reaction  temperature  in  this  circumstance  increases  with  an 
increasing  Damkohler  number  monotonically  although  an  inflection  point  is 
exhibited  in  this  relation.  The  thermochemical  intensity  of  the  reaction 
( q/e-j )  has  a  significant  effect  on  the  temperature  attained  by  the  reacting 
system. 

For  e-j  <  0.25,  the  catastrophies  of  ignition  and  extinction  are  evident. 

A  comparison  of  Figs.  6,  7a,  and  8  indicates  the  effect  of  activation  energy 
( 1/6 -j )  on  the  shape  of  the  curves.  The  ignition  Damkohler  number  is  higher  for 
higher  activation  energy  at  any  q/e-j  and  is  relatively  insensitive  to  (q/Oj ) 
at  any  activation  energy  ( l/e-j ) .  More  drastic,  however,  is  the  arrival  of 
extinction.  A  high  activation  energy  flame  is  prone  to  be  extinguished  more 
abruptly  as  evidenced  by  the  relatively  increasing  steepness  of  the  equilibrium 
branches  as  the  activation  energy  is  increased.  Additionally,  the  higher  the 
activation  energy,  the  higher  the  extinction  Damkohler  number. 

A  comparison  of  Figs.  7  a-c  shows  the  influence  of  the  fire  strength 
parameter  q  on  ignition  and  extinction  phenomena  in  a  much  greater  detail  than, 
but  quite  consistent  with.  Figs.  6  and  8.  Energetically  weaker  flames  are  more 
abrupt  in  their  response  to  perturbations  aimed  at  extinction.  Weaker  flames 
are  prone  to  be  extinguished  at  a  larger  Damkohler  number.  The  ignition  process 
however,  is  only  weakly  influenced  by  variations  in  q/e-j . 

This  last  observation  is  of  immense  importance  in  blast/fire  interaction 
studies.  For  the  case  of  medium  activation  energy  (1/e-j  =  40),  Fig.  9  shows 
the  ignition  and  extinction  Damkohler  numbers  as  dependent  upon  the  energetic 
strength  q/e-j  of  the  flames.  The  relative  insensitivity  of  ignition  is 
obvious.  More  importantly,  the  extinction  process  is  strongly  dependent  on  the 
flame's  energetic  strength. 


One  can  view  the  graph  of  Fig.  9  as  follows.  Flames  and  perturbations 
falling  in  the  area  under  the  extinction  curve  correspond  to  situations  in  which 
extinction  is  imminent;  those  falling  in  the  area  above  the  curve  are  immune  to 
extinction.  Energetically  weak  flames  (i.e.,  q  is  low)  impacted  by  strong 
blasts  (i.e.,  P  is  low),  for  instance,  fall  near  the  origin  of  the  graph  and 
are  subject  to  extinction.  Energetically  robust  flames  impacted  by  weak  blasts, 
on  the  other  hand,  fall  far  away  from  the  origin  and  are  expected  to  survive  the 
blast.  Viewed  in  another  way,  corresponding  to  any  given  fire  strength  (q/e-j) 
there  exists  a  lower  limit  blast  strength  (upper  limit  P)  below  which 
extinction  is  not  possible.  For  example,  if  q/e^  =  20,  in  P  can  not  be  lower 
than  3.5  if  extinction  is  to  be  impossible.  Alternately,  corresponding  to  any 
given  blast  strength,  there  exists  a  lower  limit  flame  strength  beyond  which 
extinction  is  impossible.  As  an  example,  for  a  blast  whose  strength  is 
described  by  tn  P  =  8,  fires  energetically  stronger  than  those  represented  by 
q/9 -j  =  7  are  capable  of  resisting  extinction. 

An  important  point  has  to  be  made  here  as  strongly  as  possible.  In  the 
above  analysis,  a  global  mechanishm  is  postulated  for  the  interaction  between  a 
flame  and  a  blast  wave.  By  mere  algebraic  manipulation  of  the  species  and 
energy  conservation  priciples,  the  concepts  of  ignition  and  extinction  are 
developed.  The  precise  gas  phase  chemical  reaction  mechanisms  are  not  always 
fully  known.  As  a  result,  global  kinetics  are  employed.  There  also  exists  a 
need  to  define  the  physical  rate  appearing  in  the  denominator  of  Damkohler 
number.  This  physical  rate  may  be  a  composite  of  both  diffusion  and  convection. 
The  exact  and  precise  relation  of  Eq.  (16)  (and  Fig.  4)  is  obtainable  only  by 
completely  solving  the  total  conservation  equations  of  mass,  momentum,  energy, 
and  species  with  appropriate  boundary  conditions.  This  solution  will  most 
probably  require  numerical  computations  due  to  the  strong  nonlinearities 
involved. 

The  simple  exposition  presented  above  is  sufficient,  to  develop  a 
methodology  to  correlate  the  existing  experimental  data,  as  shown  in  the 
following  chapter. 


Charring  Solid  Fire: 

Consider  a  piece  of  wood  burning  in  a  fire,  most  generally  in  the  company 
of  other  similar  pieces.  Combustion  of  wood  involves  at  least  two  distinct 
phases:  flaming  and  glowing. 
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The  flaming  process  is  not  unlike  that  of  liquid  fuels.  The  prime 
difference  lies  in  the  method  of  production  of  combustible  gases.  Whereas  in 
the  burning  of  a  liquid  fuel  the  fuel  vapors  are  evolved  by  a  thermodynamic 
evaporation  process  ,  they  are  produced  by  an  irreversible  thermal  chemical 
degradation  process  in  wood  burning.  While  liquid  vaporization  is  a  surface 
process  with  a  preheated  liquid  layer  under  the  vaporizing  surface,  the  solid 
pyrolysis  is  a  distributed  volumetric  chemical  kinetic  process.  In  liquid 
vaporization,  the  mass  transfer  results  in  a  regressing  surface.  In  wood 
pyrolysis,  the  exposed  surface  is  stationary  during  the  flaming  combustion 
phase.  If  excessive  gas  flow  velocities  are  involved,  erosion  and  ablation  may 
result  to  make  the  exposed  surface  regress  with  time.  ;  but  these  high 
velocities  blow  the  flame  away  to  initiate  glowing.  In  liquid  burning,  as  the 
evaporating  surface  regresses,  the  subsurface  preheated  layer  also  regresses  at 
the  same  speed  into  the  liquid,  thus  maintaining  a  temperature  field  in  the 
liquid  unaltered  in  time  in  the  frame  of  coordinates  fixed  on  the  regressing 
surface.  This  temporal  invariance  of  the  condensed  phase  heating  does  not  hold 
for  woodburning,  in  which  solid  phase  transient  conduction  continuously  and 
contiguously  tends  to  even  out  the  temperature  gradients.  Upon  continued 
flaming  combustion  of  wood,  a  carbonaceous  char  layer  is  formed  near  the 
surface.  This  layer  is  altogether  devoid  of  capacity  to  further  pyrolyze. 
Pyrolyzates  are  then  produced  within  the  solid;  they  flow  out  of  the  solid 
through  the  char  sheath,  mix  with  air  in  the  gas  phase,  combust  and  produce  the 
flame.  The  surface  char  layer  grows  in  thickness  progressively  as  the  pyrolysis 
process  penetrates  into  the  interior  of  the  solid.  Based  on  the  knowledge  of 
transient  conduction  process,  it  can  be  shown  that  the  temperature  gradients  are 
steepest  near  the  exposed  surface  and  that,  the  pyrolysis  zone  becomes 
progressively  broader  as  the  solid  is  gradually  charred  away.  If  the  solid  is 
finite  in  thickness,  a  time  will  arrive  at  which  the  entire  body  is  converted  to 
char;  no  more  gas  phase  pyrolyzates  can  be  produced;  flaming  ceases. 

The  differences  between  liquid  and  solid  fuels  undergoing  flaming 
combustion  as  described  above  are  important  in  our  attempts  to  sort  out  the 
blast/fire  interactions.  In  the  early  stages  of  wood  burning,  the  pyrolysis 
layer  is  yet  to  be  established  in  its  inward  progress  into  the  solid;  the  longer 
the  heating  and  burning  in  these  initial  stages,  the  more  sustained  the  flames 
would  become. 


As  long  as  the  flames  cover  the  burning  solid  surface,  the  oxygen 
concentration  at  the  solid  surface  is  minimal  or  zero.  If  the  flame  is  ideally 
diffusion-controlled,  oxygen  can  not  penetrate  to  the  fuel-side.  On  the  other 
hand,  if  finite-rate  gas  phase  combustion  kinetics  are  involved  or  if  the  flame 
is  mechanically  ruptured,  then  oxygen  can  reach  the  hot  solid  surface.  Glowing 
combustion  of  the  carbonaceous  solid  residue  (viz:  the  char)  is  then  expected. 

Glowing  combustion  of  carbonaceous  solids  involves  five  main  steps  in 
series:  diffusion  of  oxygen  to  the  hot  surface,  adsorption  of  this  oxygen  by 
the  surface,  reaction  of  the  adsorbed  oxygen  with  the  solid  to  form  adsorbed 
product(s),  desorption  of  the  product(s),  and  diffusion  of  the  desorbed  product 
to  the  gas  phase.  The  slowest  of  this  series  of  steps  governs  the  global 
combustion  intensity.  The  reaction  in  the  absorbed  state  can  be  expressed  as: 
C(s)  +  O2  ♦  CO2. 

Upon  closer  scrutiny,  some  investigators  propose  a  slightly  different 
mechanism,  in  which  oxygen  does  not  ever  reach  the  surface.  As  CO2  is  formed  in 
a  thin  blue  flame  in  the  gas  phase,  a  part  of  it  diffuses  to  the  solid,  gets 
absorbed,  and  participates  in  the  reduction  reaction  C(s)  +  CO2  +  2C0.  The  CO 
thus  formed  is  desorbed  to  diffuse  away  from  the  solid  and  to  encounter  the 
inward  diffusing  oxygen.  The  reaction  between  CO  and  O2  gives  a  blue  flame 
according  to  CO  +  O2  +  2CO2.  Part  of  the  CO2  thus  formed  returns  to  the  surface 
to  start  the  cycle  again.  Although  the  temperatures  of  both  the  solid  surface 
and  the  CO  flame  are  high,  the  global  carbon  consumption  rate  is  so  small  that 
the  CO  flame  is  generally  situated  quite  close  to  the  solid  surface  except  when 
the  pressure  and  the  oxygen  content  of  the  ambient  gas  are  very  low. 

The  temperature  of  charcoal  and  carbon  surfaces  is  known  to  rise  with  an 
increase  in  the  speed  with  which  air  is  blown  on  them  [9-12].  While  the  solid 
geometry,  char  grade,  air  temperature,  and  the  size  of  the  solid  body  are  found 
to  be  of  influence,  more  important  is  the  fact  that  an  increase  in  the  air  speed 
would  increase  the  surface  temperature.  This  behavior  is  in  complete  contrast 
with  what  is  expected  with  gas  phase  flaming  combustion  wherein  an  increase  in 
flow  veloicty  is  known  to  increase  thermal  dilution  as  a  consequence  of  which  the 
flame  temperature  would  fall  to  eventually  result  in  an  extinction  of  the  flame. 

In  the  case  of  heterogenous  combustion,  however,  the  very  same  mechanism 
by  which  thermal  dilution  is  realized  also  brings  about  an  enhanced  oxygen 
diffusion  to  the  reaction  site  at,  or  on,  the  surface.  It  is  this  enhanced 


oxygen  transport  which  is  reponsible  for  the  noted  dependency  of  the  reaction 
temperature  on  the  flow  velocity.  This  peculiarity  is  relevant  to  the 
blast/fire  interaction  as  described  in  the  following  simple  analysis. 


Analysis: 


Consider  a  glowing  carbon  surface.  The  oxygen  mass  fractions  in  the 
ambient  air  and  at  the  reaction  surface  are  respectively  Yq  «  and  Yqs  . 

Transfer  of  oxygen  to  the  surface  occurs  as  dependent  upon  the  flow  and  geometry 

at  a  rate  ftg"  =  ho(Yo»"YOs^  w*1ere  ®o"  15  t*ie  1,1355  02tran5|3ortecl  Per  unlt 

surface  area  per  unit  time  and  hp  is  the  mass  transfer  coefficient  (kg/m2s) 
obtainable  for  the  given  geometry  as  a  function  of  the  body  size,  diffusivity 
of  oxygen,  density  of  the  gas  phase,  and  the  Reynolds  and  Schmidt  numbers  of 
the  flow.  If  f  is  the  stoichiometric  fuel  oxygen  mass  ratio,  the  carbon  mass 
loss  rate  per  unit  area  is 


V  * f  V  W  (17) 

If  first  order  reaction  kinetics  are  assumed,  the  chemical  reaction  rate  is 
jiven  by 


ft  "  s  k  y 

"s  VOs 


(18) 


where  the  specific  reaction  rate  constant  depends  on  temperature  according  to 
Arrhenius  law. 


ks  *  k0s  e*P(-Es/RT)  (19) 

where  ES/"R  is  the  activation  temperature  for  the  glowing  reaction.  Combining 
Eqs.  (17)  and  (18), 

vo5  *  <20> 

V  ■  <21> 


ks/YhD  is  a  Damkohler  number  related  to  the  surface  glowing  reaction.  When  this 
is  very  small,  due  to  low  pressures,  temperatures  or  high  mass  transfer  rates, 
the  surface  mass  fraction  of  oxygen  departs  from  the  free  stream  value  only 
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slightly.  The  oxgen  mass  gradients  then  are  vanishingly  small.  The  process  is 

governed  by  the  slow  chemistry  so  that  &s"  =  On  the  other  hand,  if  the 

Damkohler  number  is  extremely  large  since  ks  is  very  large  compared  to  the 
relatively  sluggish  mass  transfer,  Ygs  *  Yo«./Cks/f hp]  *  0.  Severe  composition 
gradients  arise  and  the  situation  is  governed  by  the  slower  process  of  mass 
transfer  so  that  5ls"  =  fhgYo*. 

An  energy  balance  may  now  be  derived  by  combining  the  energy  and  oxygen 
species  conservation  equations  in  the  following  form  indicative  of  the 
correspondence  between  the  depletion  of  O2  and  heating  of  the  surface. 


Combining  Eqs.  (20)  and  (22),  with  Eq.  (19), 


wvT->  f<v-v 


C  pq 


where 

0  =  R  T  /E  e  s  R  T  /E 

S  S  «»  00' 


Os 


and 


This  result  is  similar  to  Eq.  (16)  which  corresponds  to  gas  phase  flaming.  In 
Figure  10,  Eq.  (23)  is  depicted  schematically.  If  the  length  of  time  the  wood 
has  been  burning  since  its  ignition  is  long,  the  char  can  be  expected  to  be  of 
high  quality  in  its  carbon  content.  Judging  from  the  experimental  data  of  Refs. 
9-12,  charcoal  glowing  temperature  is  lower  than  that  of  graphite.  As  a  result, 
one  can  expect  that  graphite  could  be  ignited  to  glow  with  lesser  effort  in 
providing  oxygen,  i.e.,  the  Damkohler  number  for  glowing  ignition  of  graphite  is 
higher  than  that  for  charcoal.  Hence  if  the  time  of  wood  burning  since  ignition 
determines  the  quality  of  char,  Eq.  (23)  for  short  and  long  burning  times  should 
appear  as  Indicated  in  Fig.  10. 

That  an  Increase  in  air  flow  velocity  (which  causes  a  decrease  in 
Damkohler  number)  will  Increase  the  glowing  surface  temperature  is  thus 
explainable.  The  main  point  is  that  the  middle  branch  of  the  solution  is 
forbidden  only  for  the  gas  phase  combustion;  for  glowing  combustion  it 
constitutes  an  allowable,  although  unstable,  solution. 


Blast/Charring  Fuel  Fire  Interaction: 


Consider  the  scenario,  based  on  one  first  enunciated  by  Stan  Martin  [13], 
involving  the  following  chain  of  events.  A  charring  solid,  such  as  wood,  is 
ignited  by  a  thermal  pulse.  The  material  burns  with  vigorous  flames  which 
progressi vley  become  stronger  as  burning  continues.  The  wood  at  the  surface 
gradually  transforms  to  char  whose  graphitic  quality  progressively  increases  as 
burning  continues  and  pyrolysis  penetrates  into  the  solid.  The  air  blast  wave 
now  arrives  to  strip,  or  blow,  the  flames  away.  The  gas  phase  adjacent  to  the 
solid  surface  becomes  rich  in  oxygen.  The  char  surface  ensues  to  glow  and  glow 
intensely  if  the  blast  winds  are  intense.  This  glowing  provides  a  strong  energy 
source  which  continuously  supplies  energy  that  is  conducted  into  the  solid  to 
cause  continued  pyrolysis.  As  the  winds  subside,  the  glowing  Damkohler  number 
increases,  the  glowing  surface  temperature  decreases  to  attain  in  time  the  frozen 
solution.  As  a  result  of  the  temperature  decrease,  the  reaction  would  further 
diminish  to  reduce  Damkohler  number  taking  the  system  down  along  the  frozen 
branch.  In  the  course  of  this  last  excursion,  if  the  system  passes  through  the 
flaming  ignition  state,  reignition  (or  reflash)  of  a  flame  would  occur.  If  not, 
the  fire  just  dies  down. 

This  scenario  is  depicted  in  Fig.  11.  The  flaming  combustion  of  a  wood 
firebrand  is  indicated  by  the  curve  marked  OlfEflf'  which  is  based  on  the  same 
ideas  as  underlying  Eq.  (16).  Flaming  ignition  occurs  from  state  If  to  state 
If1  to  achieve  a  flame  on  the  equilibrium  branch.  Branch  If Ef  is  forbidden  for 
the  flame.  As  the  firebrand  undergoes  flaming  combustion  in  a  thermochemical 
state  near  I ' p  and  as  a  blast  wave  is  inpacted  on  its  flames,  the  flame  cools 
down  towards  the  Ep  state.  If  the  blast  is  strong  enough  to  cool  the  flame  to 
Ep,  flame  extinguishment  is  imminent  in  the  same  way  as  if  the  flame  is 
supported  by  a  liquid  fuel.  There  is  one  difference,  however.  Upon 
extinguishment  of  a  liquid  fuel  flame,  the  state  E'p  is  attained.  With  charring 
fuel  flame  extinguishment,  however,  the  corresponding  post-extinction  state  is 
either  E'pi  (for  short  preburn  duration)  or  E'p2  (for  long  preburn  duration)*  on 
the  charring  unstable  middle  branch.  Upon  extinction  of  the  flame  to  arrive  at 
E'pi  or  E'p2»  continued  high  winds  may  increase  the  glowing  intensity  to  briefly 
drift  the  system  E'p2  to  A2  or  from  E'pj  to  Aj.  The  system  will  then  cool  down, 

*By  preburn  duration,  we  mean,  the  time  the  wood  has  been  burning  before  the 
blast  arrived  upon  it. 
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as  the  winds  subside,  through  the  sequence  of  states  -E'f2  “^62  or  A1  -E ' f 1 
-Iqj.  Due  to  the  corresponding  reduction  in  the  temperature  of  the  reacting 
surface,  the  chemical  kinetic  rate  itself  would  soon  begin  to  fall  from  Ig  to  0. 
In  the  course  of  this  state  transition,  the  system  may  encounter  the  flaming 
ignition  point  Ip  if  the  preburn  char  is  so  as  to  correspond  to  EqE ' p2 ^G2 J  the 
flame  would  then  reflash  to  I'p  when  the  preburn  time  is  long.  On  the  other 
hand,  if  the  preburn  time  is  short  such  as  that  corresponding  to  EqE'fiIgi,  the 
process  of  IqiO  of  cessation  of  the  reaction  would  never  pass  through  the 
flaming  ignition  state  Ip  and  hence  no  reestablishment  of  the  flame  is 
possible. 

If  the  effect  of  preburn  time  is  not  to  alter  the  graphitic  nature  of  the 
surface  char  but  to  modify  the  flammability  properties  of  the  pyrolysis  gas 
mixture,  the  flame  reflashing  may  be  explained  as  shown  in  Fig.  12.  Since  a 
longer  preburn  may  be  expected  to  improve  the  quality  of  combustibility  of  the 
pyrolyzates,  the  curve  Olf 2Ef 2^ ' f 2  represents  the  long  preburn  flaming 
situation.  Ignition  of  these  pyrolyzates  is  easier  and  the  equilibrium  flame 
temperature  is  higher.  Flaming  j>roducts  of  pyrolysis  resulting  at 

earlier  times  (i.e.,  short  prebggn  imcj)  is  representable  by  the  curve 
OlfiEfil'fi  wherein  ignition  is  relatively  difficult  and  the  heat  of  combustion 
is  relatively  smaller.  Experience  on  wood  crib  burning,  however,  leads  one  to 
believe  that  the  preburn  time  may  influence  the  quality  of  pyrolyzates  only  in  a 
minor  way. 

The  char  glowing  combustion,  invariant  with  time  of  preburn,  is  indicated 
by  curve  labelled  OIqE *fiEf2'EG.  Now  focussing  on  the  short  preburn  situation 
as  a  blast  wave  shifts  the  flame  to  Efi  extinction  occurs  to  the  char  state  E ' f x . 
As  the  winds  subside,  the  char  glowing  rate  decays  from  near  E'fi  to  Iq  and  the 
reaction  dies  from  Iq  to  0.  In  this  last  leg  of  state  oxcursion  from  Iq  to  0, 
the  system  does  not  encounter  the  flaming  ignition  state  If  1 .  Hence  this  short 
preburn  situation  involves  no  reflash  of  the  flame. 

Take  the  long  preburn  situation  in  the  same  sequence  of  arguments  and  one 
sees  that  as  the  flame  is  blown  out,  glowing  is  thwarted  by  the  subsiding  winds 
and  the  char  is  cooled  down,  the  pyrolyzate  flaming  ignition  state  If 2  is  surely 
encountered  to  result  in  a  reflash  of  the  flame  to  1 1 f 2 .  Thus  as  the  preburn 
time  Increases,  the  pyrolyzates  issued  are  more  flammable,  If  shifts  to  the  left 
from  If  1  towards  If 2 •  There  is  a  minimum  preburn  time  beyond  which  the  flaming 
Ignition  state  If  falls  on,  or  to  the  left  of,  Iq  which  corresponds  to  the 


EXPERIMENTS,  CORRELATIONS  AND  DISCUSSION 


SRI  Experiments: 

The  most  abundant  and  relevant  experiments  to  test  the  blast/fire 
interaction  hypotheses  of  the  preceding  chapter  are  those  reported  by  Martin, 
Backovsky  and  their  colleagues  at  SRI  International  [1-3].  The  experimental  set 
up  and  data  will  be  briefly  discussed  in  this  chapter  before  adapting  the  data 
for  correlations  in  the  conceptual  framework  of  ignition  and  extinction. 

The  properties  of  concern  in  the  hypothesis  may  be  conveniently  placed 
into  three  categories  respectively  pertaining  to  the  fuel,  fire,  and  the  blast. 
Among  the  fuel  properties  are:  the  pyrolysis  kinetics  of  wood  fuel  or  the 
vaporization  thermodynamics  of  liquid  fuels,  thermodynamics  of  mixing, 
thermochemistry  of  glowing  and  flaming  combustion,  oxidation  kinetics  of  glowing 
and  flaming  reactions,  and  transport  properties.  The  fire  properties  mainly 
concern  the  geometry  (shape,  size,  orientation)  of  the  fuel  bed,  the  nature  of 
flow  (forced  or  free  convection)  and  the  externally  imposed  radiant  flux.  The 
blast  properties  include  the  peak  overpressure,  time  of  positive  phase  and 
continued  radiant  pulse. 

The  SRI  experimental  facility  attempts  to  capture  these  essential 
variables  as  described  in  detail  by  Martin  and  Backovsky  in  references  1-3. 
Suffice  it  here  to  say,  the  air-driven  shocktube  is  capable  of  producing 
shockwaves  typical  of  those  produced  by  kilo-to-megaton  nuclear  explosions  in 
air.  The  peak  over  pressures  (up  to  25  psi)  and  positive  phase  durations 
(between  0.10  to  over  3.5  sec)  are  prescribable  by  the  operator  in  this 
experimental  set  up.  The  test  section  and  typical  wood  crib  are  shown  in  Figs. 
13  and  14. 

Liquid  pool  fires  are  tested  with  a  variety  of  fuels:  n-hexane, 
methanol,  kerosene,  n-pentane,  and  acetone  with  the  fuel  bed  length  ranging 
between  12  and  36  inches.  The  mean  overpressure  is  varied  between  0.91  and  7.5 
psi.  Time  of  positive  phase  duration  was  ranged  between  0.068  and  3.8  sec.  In 
some  tests,  the  preburn  time  before  the  blast  arrival  is  noted.  Barriers 
located  upstream  of  the  pool  are  also  studied  in  some  tests.  The  observations 
included  noting  whether  or  not  a  given  fire  is  extinguished  by  a  given  blast. 

Typical  wood  cribs  shown  in  Fig.  14  are  made  of  (mostly)  dry  western 
hemlock  and  (some)  redwood.  Overpressures  in  the  range  0.91  to  9.88  psi  with 


positive  phases  between  0.078-0.117  sec  are  tested.  The  cribs  are  made  of 
sticks  either  3/8"  or  3/4"  thick.  A  tray  of  propyl  alcohol  under  the  crib  is 
used  to  prompt  the  crib  ignition  following  which  the  crib  is  allowed  to  burn  for 
anywhere  from  about  60  -  180  sec  before  the  blast  wave  is  imposed  on  the  fire. 

There  are  also  reported  eleven  tests  in  which  trays  of  shredded  blotter 
or  filter  paper  act  as  fuel  beds. 

Correlation: 

We  first  acknowledge  that  it  is  always  a  precarious  task  to  undertake 
interpretation  of  other  investigators'  experimental  data.  Only  the  experimentor 
knows  all  the  subtle  details  of  the  set-up,  conduct  and  behavior  of  an 
experimental  run.  However,  opportunity  exists  to  improve  the  current  knowledge 
of  blast/fire  interaction  by  testing  the  hypothesis  of  this  report.  We  take  the 
SRI  data  at  their  face  value  as  documented  in  references  1-3  and  as  understood  in 
the  sporatic  discussions  between  the  experimentors  and  the  present  author. 

Having  made  note  of  this  limitation,  we  now  seek  to  gather  the  information 
required  to  adapt  the  experimental  data  to  the  hypothesis  of  the  preceeding 
chapter . 

Class  B  Fires: 

Our  intent  is  to  present  the  SRI  data  in  the  format  of  Fig.  9.  In  order 
to  transform  the  coordinates  into  nondimensional  form,  the  following  details  are 
to  be  noted. 

q  =  R  hcY/\i/E  Cng.  The  universal  gas  constant  requires  no  comment. 
Enthalpy  of  combustion  for  various  liquid  fuels  is  well  known.  The  mass  fraction 

at  the  fuel  boundary  depends  on  the  pressure  as  well  as  temperature  at  the 

boundary.  Prior  to  the  blast  arrival  the  surface  concentration  will  correspond 

to  equilibrium  vaporization.  Immediately  following  the  blast  arrival,  the  flame 

will  be  disturbed;  there  will  be  significant  mixing.  The  rise  in  temperature  due 

to  shock  compression  will  be  too  rapid  for  the  condensed  phase  to  promptly 

respond.  We  take  here  the  fuel  concentration  to  be  same  as  that  in  stable 

equilibrium  diffusional  combustion.  For  most  hydrocarbon  fuel  pools  burning  in 

air,  the  fuel  concentration  at  the  surface  can  be  estimated  to  be  about  1  kg/m^. 

This  quantity  is  obviously  a  function  of  the  saturation  vaporization  process  in 

which  such  parameters  as  the  liquid  temperature,  latent  heat  of  vaporization  play 
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a  role.  The  precise  delineation  requires  a  fundamental  study  of  liquid 
vaporization,  such  a  study  being  out  of  the  scope  of  the  current  project.  The 
quantity  Y/^  thus  would  be  taken  as  l(kg/m^)/p (kg/m^)  where  p  is  the  air  density 
at  a  pressure  representative  of  the  conditions  immediately  following  the  blast 
arrival.  The  specific  heat  of  gases  will  be  taken  to  be  a  constant  equal  to  1666 
J/kgK  irrespective  of  the  temperature  and  composition.  This  again  is  an 
approximation  without  which  the  application  becomes  horrendous.  Literature  [7] 
on  combustion  of  hydrocarbons  in  air  indicates  that  the  global  activation  energy 
E  ranges  between  about  150  to  250  MJ/kgmole  depending  upon  the  fuel  an  the 
combustion  situation.  In  studies  not  seeking  to  elicit  the  actual  elementary 
chemical  reaction  steps  culminating  in  the  global  behavior,  it  is  fairly  safe  to 
take  E  to  be  approximately  200  MJ/kgmole,  independent  of  the  hydrocarbron .  One 
thus  estimates  q  as 

q  «  2.8x10-7  hc(in  J/kg)/P(in  psia)  (24) 

The  units  are  inelegantly  mixed  because  hc  is  generally  known  in  SI  units  and  the 
experimental  tabulations  give  pressure  in  psi . 

P  =  kQt/u.  We  take  the  velocity  u  to  be  the  particle  velocity  produced 
by  the  blast.  By  the  relation  given  on  page  6,  with  y  =  1.4  for  air  whose 
pressure  prior  to  shock  arrival  is  Po  =  14.7  psia,  and  speed  of  sound  in  ambient 
air  co  =  1076  ft/s,  we  find  u(in  ft/sec)  =  1076  p/[20.58x(l+0.0583p)l/2]  where  p 
is  the  peak  overpressure  in  psi.  Furthermore,  the  global  collision  frequency  kg 
of  the  assumed  first  order  oxidation  reaction  is  taken  to  be  5xlolO(l/sec) 
independent  of  the  fuel  within  the  same  restraints  of  approximation  as  the  global 
activation  energy  discussed  above.  One  would  thus  obtain 

P  *  1.5x10^(1+0.0583  p  (in  psi ) 1/2  *(in  inches)/p  (in  psi)  (25) 

again  the  units  being  written  for  the  convenience  of  adoption  to  experimental 
data,  l  is  the  fuel  bed  length. 

Tables  1-4  are  directly  taken  from  references  1-3  to  report  the  class  B 

fire  data  involving  no  barriers.  Forty  six  different  tests  are  involved  in  these 

tables  on  five  different  fuels  with  three  tray  lengths.  The  blast  strength 

parameter  P  and  fire  strength  parameter  q  are  estimated  for  each  rur  according 

to  the  preceding  explanation  and  columns  are  added  to  the  original  tables.  The 
data  are  shown  plotted  in  Fig.  15  in  coordinates  suggested  by  our  hypothesis  of 


Fig.  9  delineating  the  regimes  in  which  the  flame  would  be  extinguished  from 
those  in  which  it  would  sustain  the  blast.  The  closed  symbols  correspond  to 
flames  which  are  extinguished  while  the  open  ones  stand  for  flames  which  resisted 
extinction. 

The  boundary  curve  shown  in  Fig.  15  is  precisely  the  same  as  the  one  in 
Fig.  9  except  that  it  is  shifted  upwards  by  a  factor  of  10^.  This  is 
permissible  since  our  entire  theoretical  effort  is  not  meant  to  be  a  development 
of  a  prediction  capability.  Instead  its  goal  is  to  develop  concepts  and  general 
forms  of  the  nondimensional  descriptions  of  the  paramters.  Thus  shifting  the 
curve  of  Fig.  9  so  that  it  passes  through  q  =  0.55  and  P  =  25,  the  boundary 
between  the  extinction  domain  and  nonextinction  domains  is  obtained. 

With  the  exception  of  the  four  anomalous  experiments  in  which  extinction 
occurred  when  it  is  not  expected*,  the  hypothesis  bears  well  in  separating  the 
conditions  suitable  for  extinction  from  those  unsuitable.  In  fact  the  goodness 
of  this  separation  is  especially  impressive  if  one  recalls  that  the  hypothesis 
is  based  on  a  theoretical  model  of  remarkable  simplicity. 

The  global  combustion  chemistry  kinetics  are  not  well-known.  Transport 
properties  are  poorly  defined.  The  mixing,  cooling,  and  dilution  following  the 
arrival  of  the  blast  wave  is  certain  to  be  hopelessly  complicated  for  any  level 
of  theoretical  description.  There  will  undoubtedly  be  alterations  in  the 
vaporization  thermodynamics  as  the  blast  effects  of  imposed.  The  surface  fuel 
vapor  concentration  assumption  is  admittedly  an  over-simplification.  Not 
withstanding  all  the  rough  spots  in  the  evolution  of  the  theoretical  hypothesis 
underlying  the  map  of  Figs.  9  and  15,  the  success  is  indeed  gratifying. 

Attention  is  drawn  to  the  fact  that,  in  preparing  Fig.  15,  the  length  of 
positive  phase  duration  is  not  at  all  invoked.  That  even  this  omission  did  not 
influence  the  validity  of  the  hypothesis  underlying  Fig.  9  indicates  that  if  the 
fire  did  not  get  seriously  disturbed  towards  extinction  in  the  very  first  few 
milliseconds  of  blast  impact,  it  most  probably  is  going  to  survive  the  blast. 

This  is  an  interesting  inference,  for  it  suggests  that  the  future  predictive 
theories  may  perhaps  be  based  on  a  description  of  the  events  promrtly  surrounding 
the  very  initial  stages  of  the  blast/fire  interaction. 

The  SRI  data  also  include  some  runs  in  which  a  1.75"  barricade  wall  is 
placed  upstream  of  the  fire.  Although  theoretical  solution  of  the  problem 

♦These  anomalous  runs  are  run  #7  of  Table  1  on  hexane,  run  #2  of  Table  3  on 
kerosene  and  runs  #8  and  10  of  Table  3  on  pentane. 


describing  the  influence  of  barricades  on  blast/fire  interaction  requires 
detailed  analysis  and  calculations,  one  can  expect  barricades  to  generally 
stabilize  fires.  As  the  flow  occurs  over  a  barricade  step,  strong  recirculation 
of  flow  is  expected  behind  it.  Such  a  recirculation  will  in  effect  slow  the  flow 
down  (or  increase  the  flow  path  length  and  flow  time).  The  net  result  is  that 
the  presence  of  a  barricade  will  increase  the  nonextinction  domain  of  the  P-q 
map  by  displacing  the  boundary  curve  downwards. 

The  SRI  barricaded  pool  fire  data  from  Reference  2  is  reported  in  Tables  5 
and  6  with  columns  added  for  P  and  q.  The  data  are  plotted  in  Fig.  16.  Except 
for  the  anomalous  run  #64  of  Table  6  on  hexane,  the  downward  shift  of  the 
boundary  curve  is  evident  both  for  hexane  and  methanol  fires.  Note,  however  that 
the  height  of  barricades  tested  (1.75")  is  relatively  small  compared  to  the 
length  of  even  the  shortest  fuel  bed.  Such  short  barricades  at  best  are  expected 
to  modify  the  flow  profiles  without  producing  significant  stable  recirculation. 
Experiments  in  which  the  barricade  height  is  comparable  in  magnitude  to  the  pool 
length  appear  to  be  desirable  to  simulate  the  influence  of  a  window-sill  in  the 
interaction  of  a  blast  with  a  fire  in  a  room.  No  systematic  extinction- 
nonextinction  boundary  curves  can  now  be  drawn  on  the  basis  of  the  available 
short  barricade  data  presented  in  Fig.  16. 

In  Reference  3  Backovsky  draws  the  conclusion,  based  on  high-speed 
photographic  observations,  that  the  most  plausible  mechanism  of  class  B  fire 
extinction  by  blast  is  one  based  on  displacement  of  the  flame  from  the  fuel 
surface.  This  is  not  inconsistent  with  the  hypothesis  of  this  report  in  which 
the  definition  of  extinction  evolved  to  be  an  annihilation  of  the  flame  by 
whatever  specific  details.  Excessive  cooling  and  dilution  to  levels  prohibiting 
flame  reflashing  from  downstream  is  within  the  realm  of  our  hypothesis.  Through 
the  success  apparent  in  Figs.  15  and  16,  it  is  fair  to  say  that  expectable  trends 
of  extinction  behavior  are  achieved  as  dependent  upon  the  fuel  type,  pool  size, 
and  blast  strength,  with  as  well  as  without  the  barricades. 

Class  A  Fires: 

The  SRI  experiments  on  blast  interaction  with  wood  crib  fires  indicate 

that  the  length  of  time  a  crib  is  permitted  to  burn  before  the  blast  imposition 

(i.e.,  the  'preburn  time')  has  a  strong  effect  on  the  extinction  process.  Cribs 

with  longer  preburn  time  require  stronger  blasts  for  extinction.  Even  more 
interestingly,  the  experiments  show  that  if  the  preburn  time  exceeds  a  certain 


minimum  value  (-170s  corresponding  to  a  minimum  fractional  weight  loss  »  35%), 
the  fire  sustains  the  impact  of  any  blast  (in  the  tested  1-10  psi  overpressure 
range)  irrespective  of  the  blast  strength. 

In  attempts  to  represent  the  blast/fire  interaction  for  class  A  fires  on 
the  same  basis  as  class  B  fires  presented  above,  several  questions  arise.  These 
questions  mainly  pertain  to  the  pyrolysis  of  wood  and  wood-like  solids  as  they 
burn  in  the  format  of  a  crib.  Specifically  how  does  the  preburn  time  make  the 
fire  progressively  stronger  and  eventually  strong  enough  to  become  immune  to 
extinction?  Said  another  way,  how  does  one  nondimensionl ize  the  preburn  time  to 
obtain  a  parameter  which  is  capable  of  serving  as  a  measure  of  the  fire's 
energetic  strength?  One  can  surmise  that  continued  heating  of  wood  will 
gradually  enhance  the  combustibility  quality  of  the  pyrolyzates.  But  this 
speculation  is  dismissed  by  the  observation  that  the  crib  burning  rate  attains  a 
steady  maximum  value  in  a  reasonably  short  time.  Indeed,  it  appears  that  no 
argument  based  on  gas  phase  combustion  phenomena  yields  a  viable  basis  to  explain 
the  SRI  experimental  trends. 

Therefore,  we  turn  to  the  solid  phase  phenomena  of  transient  conduction 
and  pyrolysis.  As  crib  burning  continues  during  the  preburn  time,  the  solid  is 
continuously  heated  up  due  to  conduction.  Also  the  pyrolysis  process  develops 
within  the  solid  transiently  to  evolve  into  an  inward  propagating  wave.  Our 
following  treatment  is  based  on  the  premise  that  the  transient  heating  process 
improves  the  blast  resistance  of  a  crib  fire  as  progressively  more  energy  is 
stored  within  the  fuel  element  to  sustain  a  certain  inertia  to  the  pyrolysis  and 
combustion  process.  An  appropriate  time  scale  for  this  effect  appears  to  be  the 
burning  time  constant  t®  which  represents  the  time  taken  by  a  fuel  element  to 
completely  char.  The  nature  of  t^  is  discussed  in  later  paragraphs. 

If  the  ratio  of  preburn  time  to  the  burning  time  constant  represents  the 
gradually  increasing  blast  resistance  of  the  crib  fire  by  the  above  argument, 
what  factors  determine  the  critical  preburn  time  beyond  which  the  fire  becomes 
blast-proof?  The  answer  to  this  question  is  sought  under  the  premise  that  at 
this  critical  time,  the  surface  of  the  fuel  element  is  completely  charred  so  that 
at  times  longer  than  this  the  pyrolosis  wave  is  fully  developed  to  propagate  into 
the  solid  thus  gradually  increasing  the  thickness  of  the  surface  char  layer. 

Based  on  pyrolysis  kinetic  data,  the  critical  preburn  time  can  be  estimated. 


Burning  Time  Constant  t®: 


Backovsky  et  al  [2,3]  give  an  excellent  summary  of  the  state-of-the-art 
knowledge  about  the  burning  behavior  of  wood  cribs.  Most  notable  in  this  summary 
are  the  facts  that:  (a)  the  crib  burning  quickly  attains  a  steady  state  in  which 
more  than  the  middle  50%  of  the  total  mass  loss  would  occur  such  that  mass 
decreases  linearly  with  time;  (b)  this  steady  state  burning  rate  depends  upon 
both  the  stick  thickness  and  the  packing  density  of  the  crib;  (c)  for  loosely 

packed  cribs  the  steady  burning  rate  is  given  by  m  =  K  A  //E~  where  m  is  the 
rate,  As  is  the  total  surface  area  of  the  exposed  wood  in  the  crib  and  b  is  the 
stick  thickness.  The  constant  K  depdnds  upon  the  type  and  moisture  content  of 
the  wood.  K  -  1.05x10-3  gm/cml*3  sec  for  pine  cribs  with  moisture  content  in  the 
range  5-10%;  and  (d)  the  feedback  from  the  flame  in  the  plume  over  the  crib 
(i.e.,  'topflame')  is  not  a  significant  factor  in  determining  ,>e  steady  rate  of 
crib  burning. 

The  square-root  dependence  of  burning  rate  on  stick  thickness  is  empirical. 
Kanury  [14]  theoretically  deduced  the  time  to  completely  char  a  cylindrical  wood 
element  under  well -ventilated  crib  conditions  as 

tO  *  A  b  +  B  b2 

where  b  is  the  element  thickness  and  A  and  B  are  constants  defined  as  below. 

A  =  [L(po-pc)  +  PoCpo(Tp'-To)/2]/2q" 

B  =  1/(32  a0) 

where  L  is  the  latent  heat  of  pyrolysis  of  wood  (*  300  J/g  volatile  production), 

PO  is  original  wood  density  (»  0.6  g/cm3),  pc  is  final  char  density  (  »  0.15 

g/cm3),  CpO  is  original  wood  specific  heat  (»  2J/g  K),  Tp *  is  temperature  at  which 

wood  pyrolysis  is  significant  (»  600  K),Tq  is  initial  temperature  («  300  K),  ao  is 

the  wood  thermal  diffusivity  (»  1x10"°  cnr/s)  and  q"  is  the  heat  flux  involved  in 

the  mutual  interaction  among  elements  within  the  crib  (•  0. 5-2.0  J/cm2s).  Noting 
that  t®  is  related  to  steady  burning  rate  iti  through  the  crib  initial  mass  m^ 

according  to  t^  -  mn/m.  The  empirical  relation  highlighted  by  Backovsky  leads  to 


t°  «  bl*5.  However,  our  theoretical  relation  indicates  that  while  thin  fuel 
elements  have  a  burning  time  t°  proportional  to  b,  thick  elements  are  consumed  as 
t®  «  b2.  While  the  heat  flux  and  pyrolysis  energetics  (and  kinetics)  play  a 
dominant  role  in  thin  element  pyrolysis,  transient  conduction  within  the  solid 
predominantly  governs  the  thick  element  pyrolysis.  The  effect  of  moisture  in  the 
solid  can  be  taken  into  account  by  increasing  the  heat  of  pyrolysis  L. 

p 

For  reasonably  dry  wood  then,  taking  q"  to  be  0.6  J/cm  s, 

t°  -  262  b  +  31  b2 

where  t°  is  in  seconds  and  b  is  in  centimeters.  The  fire  strength  q  then  is 
proportional  to  the  ratio  of  preburn  time  to  t®.  A  plot  of  the  blast  strength 
paramter  P  against  this  ratio  is  expected  to  differentiate  conditions  of  crib 
flame  extinguishment  from  those  of  nonextinguishment.  By  the  experience  obtained 
from  Figures  16  and  17,  since  the  boundary  curve  appears  to  be  same  for  different 
fuels  and  since  the  wood  pyrolyzates  are  mostly  a  mixture  of  several  hydrocarbons 
and  carbohydrates,  there  is  reason  to  believe  that  the  boundary  curve  for  wood 
flame  extinction  is  approximately  same  as  that  for  class  B  fires.  As  a  result, 
q  for  charring  fuel  fires  is  denoted  by  t  and  defined  as  t  =  2.7x  preburn 
time/t0.  The  factor  2.7  is  obtained  from  trial  and  error  to  account  for  possible 
differences  between  class  A  and  class  B  fires.  A  justification  for  this  is 
yet  to  be  developed. 

Surface  Charring  Kinetic  Time: 

Wood  pyrolysis  is  generally  known  to  be  an  extremely  complex  chemical 
kinetic  process.  However,  it  is  customary  in  most  theoretical  analyses  to  take 
the  pyrolysis  process  to  be  of  order  unity  with  the  preexponential  factor  and 
activation  energy  respectively  in  the  ranges  10^-10^  1/s  and  100-200  kJ/gmole. 
Taking  them  to  be  10^  1/s  and  140  kJ/  gmole,  we  can  estimate  the  time  taken  by 
the  surface  to  completely  pyrolyze  to  be  about  1.6x10^  s  at  T  =  500  K,  464  s  at 
600  K,  7.2  s  at  700  K.  The  surface  temperature  of  transiently  heated  wood  rises 
nearly  as  the  square  root  of  heating  time.  Pyrolysis  experiments  reported  by 
Kanury  [15]  Indicate  that  the  surface  takes  about  2  to  3  minutes  of  exposure  to 
completely  transform  to  char  in  a  fire  situation.  Thus,  the  above  estimation 
leads  to  a  surface  temperature  of  burning  wood  to  be  about  620  K  at  which  the 
time  to  completely  pyrolyze  is  about  170  s. 
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The  critical  preburn  time  beyond  which  the  crib  fire  becomes  blast-proof 
is  noted  in  the  SRI  experiments  to  be  also  of  the  order  of  170  s.  Thus  there 
seems  to  be  reason  to  believe  that  this  critical  preburn  time  corresponds  to  the 
complete  charring  of  the  wood  surface.  Upon  producing  a  char  sheath  with 
preburning  longer  than  this  critical  time,  the  pyrolysis  process  is  confined  to 
the  interior  of  the  solid  to  become  relatively  insensitive  to  the  blast  effects 
in  its  immediate  response.  Additionally,  the  surface  char  has  then  become  a  more 
effectively  glowing  fuel  to  respond  to  the  blast  effects  immediately  in  producing 
conditions  of  flame  reflashing  by  the  hypothesis  postulated  earlier.  One  would 
thus  expect  the  flaming  boundary  curve  to  be  transformed  to  a  vertical  line 
abruptly  at  t  =  2.7x  critical  preburn  time/t^.  Since  both  the  critical  preburn 
time  and  t°  depend  upon  the  crib  stick  thickness,  this  critical  t  may  depend  upon 
the  stick  thickness.  Critical  t  is  estimated  to  be  0.85  and  0.65  respectively 
for  3/8  inch  and  3/4  inch  stick  cribs.  However,  since  scaling  should  alleviate 
this  dependency  on  stick  thickness  and  since  the  relation  between  q  and  t  is  yet 
to  be  justified  by  a  reasoned  hypothesis,  we  take  the  critical  t  to  be  0.75, 
average  of  the  values  for  the  two  thicknesses.  Even  if  tentatively,  this 
critical  t  is  presumed  to  be  independent  of  stick  thickness. 

Tables  7-10  are  directly  taken  from  References  2  and  3  to  report  class  A 
fire  data  with  columns  added  for  the  blast  strenth  parameter  P  and  the  fire 
strength  parameter  t.  Thirty  nine  tests  with  3/4-inch  stick  cribs,  thirty  four 
tests  with  3/8-inch  sticks  and  ten  tests  with  trays  of  shredded  paper  are 
involved  in  this  series.  Three  3/4-inch  redwood  stick  cribs  are  also  tested  and 
data  given  in  Table  4.  In  nine  of  these  86  tests  flame  reflashing  is  noticed  as 
the  blast  subsides.  Although  glowing  of  different  intensities  is  observed  in 
over  a  dozen  tests,  not  all  glowing  tests  experienced  the  flame  reflashing.  In 
some  tests,  the  crib  is  noted  to  have  burned  poorly.  In  some  others,  the  shock 
tube  diaphragm  ruptures  poorly.  In  yet  others,  the  reflashing  occured  as  the 
test  section  is  opened  up  to  remove  the  extinguished  crib  from  the  tube.  In  two 
of  the  shredded  paper  tray  tests,  reflashing  is  found  to  occur  when  air  is  fanned 
over  the  glowing/smoldering  fuel  bed. 

Suffice  It  to  note,  unlike  class  B  fires,  the  crib  fire  experiments  are 
extremely  difficult  to  control.  Although  a  great  deal  of  care  is  excercised  by 
the  SRI  team  In  the  conduct  of  these  experiments  and  although  much  valuable 
information  is  developed,  there  persist  numerous  ambiguities  some  of  which  may  be 
Impossible  to  be  ever  resolved.  Keeping  this  in  mind,  the  data  of  Tables  7-10 
are  shown  plotted  in  Fig.  17. 
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Although  the  correlation  of  Fig.  17  appears  to  be  less  clear  cut  than  the 
Class  B  fire  correlations  of  Figs.  15  and  16,  the  success  of  the  hypothesis  is 
evident.  All  the  shredded  paper  tray  fires,  one  3/4-inch  redwood  crib  fire,  and 
eight  3/8-inch  hemlock  stick  crib  fires  defy  the  hypothesis  of  this  report  by 
failing  to  resist  the  blast  extinction.  Similarly,  ten  (3/4-inch  hemlock  crib) 
tests  fail  to  extinguish  where  they  are  expected  to.  The  most  probably  reason 


underlying  these  delinquent  nonextinctions  is  that  the  thicker  sticks  in  the  crib 
provide  more  effective  recirculation  pockets  conducive  for  stabilizing  the 
flames.  Nonetheless,  it  can  be  seen  that  the  extinctions  lie  mostly  within  the 
extinction  envelope  and  the  nonextinctions,  without  it.  Most  of  the  reflashing 
is  noticed  at  t  larger  than  0.75  which  may  be  recalled  to  correspond  to  complete 
charring  of  the  solid  surface.  Almost  all  of  these  reflashed  tests  involve  3/8- 
inch  stick  cribs. 

The  shredded  paper  experiments  seem  in  their  majority  to  defy  the 
hypothesis  of  this  report.  This  behavior  is  perhaps  understandable  on  the  basis 
that  smoldering  of  the  shredded  fuel  bed  is  in  many  respects  different  from  the 
burning  of  well -ventilated  (loosley  packed)  stick  cribs.  Furthermore,  the 
mechanical  integrity  of  the  shredded  fuel  beds  is  expected  to  be  so  poor  as  to 
enhance  the  cooling  of  the  fuel  elements  by  the  blast  winds.  Flame  holding 
recirculation  zones  are  also  minimal  in  these  beds.  Based  on  such  factors,  it 
appears  reasonable  to  expect  shredded  fuel  beds  to  be  more  susceptible  to 
extinction  by  blast. 

No  particular  pattern  is  discernible  from  the  present  correlation  as  to 
the  reflashing  mechanism.  Thus  the  hypothesis  postulated  earlier  can  not  be 
systematically  tested  at  present. 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  Damkohler  number  P  and  the  fire  strength  parameter  q  appear  to 
constitute  the  major  parameters  in  describing  the  mechanism  of  'blowing-away'  of 
flames  due  to  the  impact  of  a  blast  wave.  While  the  Damkohler  number 
represents  inverse  of  blast  stength,  the  fire  strength  parameter  represents 
mainly  the  heat  of  combustion.  Based  on  an  algebraic  model,  a  concept  is 
developed  to  identify  a  map  of  P  versus  q  on  which  the  domain  of  conditions 
conducive  to  fire  extinction  can  be  separated  from  conditions  unconducive.  The 
hypothesis  is  tested  by  the  SRI  shocktube  experimental  data  on  a  variety  of  class 
B  and  class  A  fires.  Flames,  on  both  liquid  and  solid  fuels,  appear  to  be  blown 
out  by  the  same  mechanism.  As  shown  in  Figs.  15  and  17,  the  P  required  to 
achieve  extinction  is  lower  for  flames  of  higher  q. 

The  influence  of  an  upstream  barricade  is  shown  to  render  the  liquid  fuel 
flames  more  blast-resistant  so  that  the  nonextinction  region  of  the  P(q)  map  is 
increased.  No  distinct  boundaries,  however,  can  be  drawn  now  for  the  barricade 
effect.  Experiments  are  required  in  which  taller  barricades  are  systematically 
tested.  Theory  describing  the  mechanism  of  flame  stabilization  by  recirculation 
has  to  be  developed  to  compare  with  these  experiments. 

The  enthalpy  of  combustion  of  the  fuel,  fuel  bed  length,  barricades  and 
shock  wave  strength  are  the  most  important  variables  in  the  blast/class  B  fire 
interaction  process.  The  fuel  vaporization  thermodynamics,  including  enthalpy  of 
vaporization,  the  nature  of  wick  holding  the  fuel  in  the  pool,  and  combustion 
kinetic  parameters  appear  to  have  little  or  no  effect  on  the  response  of  class  B 
fires  to  blasts.  This  is  a  significant  point  since  such  properties  as  those 
pertaining  to  global  oxidation  kinetics  can  not  be  expected  ever  to  be  known  in 
reliable  and  quantitative  detail. 

The  hypothesis  of  blast/fire  interaction  postulated  and  tested  with  pool 

fire  data  in  this  report  is  capable  of  describing  fire  extinguishment  viewed  in 

different  ways.  For  instance,  if  extinguishment  is  defined  as  annihilation  of 
the  flame,  it  includes  both  the  degeneration  of  a  stationary  flame  and 
'nonreappearance'  of  a  displaced  flame.  The  flame  displacement  mechanism 
advocated  by  the  SRI  group  thus  is  completely  consistent  with  the  present 
hypothesis.  The  success  evident  in  the  correlations  of  Fig.  15  is  probably  the 
proof  of  this  consistency. 

Having  correlated  the  12  to  36  inch  long  liquid  pool  fires,  we  believe  the 

next  task  in  studying  class  B  fires  is  to  develop  ways  of  scaling  the  fires  up  to 
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several  tens  or  hundreds  of  feet  in  size.  Development  of  these  scaling  rules 
requires  consideration  of  such  additional  factors  as:  thermal  radiation 
feed-back  from  the  larger  flames  to  the  fuel  surface,  effect  of  the  size  of  the 
flames  on  the  uniformity  with  which  the  blast  would  perturb  them,  accentuation  of 
the  recirculation  phenomena,  and  others.  Based  on  the  definition  of  P  ,  the 
critical  fuel  bed  length  is  nearly  proportional  to  the  peak  overpressure  and  the 
particle  velocity.  (See  Eq.  (25)).  This  conclusion  is  similar  to  that  reached 
earlier  by  Backovsky  et  al  in  Ref.  2. 

Turning  now  to  class  A  fires  as  represented  by  cribs  of  wood  sticks,  we 
note  that  the  mechanism  of  flame  extinction  appears  to  be  the  same  as  that  of 
class  B  fires.  The  extinguished  state,  however,  involves  charcoal  embers  capable 
of  undergoing  a  metastable  glowing  combustion  process.  This  glowing  combustion 
process  supplies  sufficient  energy  to  the  solid  to  continue  production  of 
pyrolyzates.  As  the  blast  effects  subside,  the  cooling  (and  yet  pyrolyzing) 
embers  may  pass  through  a  thermochemical  state  corresponding  to  reflashing  of 
flame  in  the  pyrolyzate  and  air  mixture. 

Based  on  the  knowledge  of  the  composition  of  pyrolysis  products,  and  of 
the  relative  independence  of  the  extinction  boundary  of  class  B  fires  on  the 
specific  nature  of  the  fuel,  there  seems  to  be  no  compelling  reason  to  believe 
that  the  flame  extinction  boundary  for  crib  fires  will  be  different  from  that  for 
class  B  fires.  The  correlation  of  P  versus  q  to  confirm  this  idea  is  given  in 
Fig.  17.  This  correlation  is  not  as  clear-cut  as  that  of  liquid  fuel  flames. 

The  reasons  for  this  can  be  any  combination  of  such  complications  as 
following,  (a)  Although  it  appears  logical  to  expect  the  preburn  time  to  be  an 
indirect  measure  of  the  fire  strength  q,  the  exact  connection  is  not  clear  at 
present.  Further  study  is  required  of  this  issue  to  establish  the  role  of 
transient  conduction  (with  pyrolysis)  in  rendering  persistence  and  strength  to  a 
crib  fire,  (b)  In  a  crib,  opportunity  exists  for  the  flames  to  find  intense 
recirculation  zones  behind  individual  sticks.  These  recirculation  zones  are  more 
effective  if  the  sticks  are  thicker.  Because  recirculation  tends  to  stablize  the 
flames  and  thus  to  shift  the  extinction  boundary  downwards,  it  is  reasonable  to 
expect  nonextinctions  where  extinction  would  be  expected  in  the  absence  of 
recirculation.  That  this  is  correct  is  visible  in  Fig.  17,  especially  in 
association  with  the  (3/4  inch)  thick  stick  fires.  The  future  study  recommended 
above,  of  flame  stablization  by  recirculation,  will  aid  in  understanding  the 
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effect  of  stick  thickness  on  extinction.  The  high  surface  temperature  of  the  wood 
sticks  in  a  crib  will  be  an  important  additional  parameter  to  be  considered  in 
this  connection,  (c)  The  surfaces  at  different  locations  in  the  crib  would  char 
nonuniformly  in  all  crib  fires.  The  distribution  of  glowing  patches  which  may 
expediently  pilot  the  reignition  of  blown-out  flames  is  thus  a  serious  ambiguity 
of  all  crib  fires.  Resolution  of  this  ambiguity  requires  further  work,  (d) 
Compared  to  liquid  pools,  the  crib  fuel  beds  are  volumetric  by  two  different 
counts,  first  by  the  crib  structure  itself  and  second  by  the  conduction-controlled 
interior  of  the  individual  sticks.  While  the  flames  over  pools  are  all-important 
in  determining  the  pool  burning  rate,  the  flame  over  the  crib  (i.e.,  the 
'top-flame')  is  of  little  or  no  significance  compared  to  the  flames  within  the 
crib,  (e)  Thermal  radiation  from  the  flames  to  the  fuel  surface  is  more  or  less 
unimportant  in  the  class  B  fires  in  the  size  range  tested.  In  contrast,  thermal 
radiation  plays  a  dominant  role  in  the  life,  strength  and  death  of  a  crib  fire. 
Additionaly  study  is  required  to  gain  a  better  understanding  of  this  role. 

The  data  correlations  of  Fig.  17  also  point  out  the  possibility  that  the 
basic  burning  process  of  shredded  paper  fuel  beds  may  be  altogether  different  from 
that  of  cribs.  It  may  be  closer  to  that  of  liquid  fuels,  but  this  is  yet  to  be 
examined  by  additional  study. 

The  flame  reflashing  mechanism  postulated  in  this  report  can  not  at  present 
be  tested  by  a  correlation  of  the  sort  shown  in  Fig.  17.  It  is  surprising, 
however,  to  note  how  few  experiments  in  the  SRI  collection  in  fact  show  a  distinct 
reflashing.  In  as  much  of  charcoal  glowing  is  expected  to  play  an  important  role 
in  this  reflashing  process,  it  appears  important  that  a  detailed  study  be 
undertaken  of  the  effect  of  blast  on  glowing  wood  surfaces.  The  nature  of  glowing 
as  it  is  distrubed  by  the  blast  variables,  the  implications  in  continued 
pyrolysis,  and  possible  transition  of  glowing  to  flaming,  are  not  only  important 
in  understanding  the  reflashing  process  but  also  exciting  from  the  view-point  of 
science. 

Turning  to  scaling,  not  withstanding  the  complicated  nature  of  wood  crib  fire 
burning  discussed  above,  crib  fires  are  probably  easier  to  be  scaled  up  since 
their  burning  is  more  or  less  locally  controlled.  Oxygen  depletion  in  the  depths 
of  the  crib  may,  however,  become  a  crucial  factor  in  such  scaling  endeavors. 
Thermal  radiation  effects,  of  course,  are  critically  required.  One  important 
conclusion  of  the  present  study  is  that  wood  crib  burning  with  and  without  blast 
interaction  requires  much  further  study  before  results  useful  in  efforts  of 
mitigation  can  be  obtained. 


REFERENCES 


1.  Martin,  S.B.,  "Experiments  on  Extinction  of  Fires  by  Airblast,"  Annual  Report, 
Contract  DCPA01-79-C-0245,  FEMA  WU  2564A,  SRI  International,  Menlo  Park,  CA 
(January  1980). 

2.  Backovsky,  J.,  Martin,  S.,  and  McKee,  R.,  "Blast  Effects  on  Fires,  "Annual 
Report,  Contract  0CPA01-79-C-0245,  FEMA  WU  2564A,  SRI  International,  Menlo 
Park,  CA  (December  1980). 

3.  Backovsky,  J.,  Martin,  S.B.,  and  McKee,  R.,  "Experimental  Extinguishment  of 
Fire  by  Blast,"  Final  Report,  Contract  EMW-C-0559,  FEMA  WU  2564A,  SRI  Inter¬ 
national,  Menlo  Park,  CA  (March  1982). 

4.  Glasstone,  S.,  and  Dolan,  P.,  The  Effects  of  Nuclear  Weapons,  U.S.  Departments 
of  Defense  and  Energy,  USGPO  stock  No.  008-046-00093-0,  (1977). 

5.  Rockett,  J.A.,  "Outline  of  Fire  Development  Models  and  Their  Status,"  pp.  D-l 
to  D-4,  Proceedings  of  Blast/Fire  Interactions,  Asilomar  1979  Conference.  R.S. 
Alger,  S.B.  Martin,  Editors,  SRI  International  Report  for  DCPA,  Contract  No. 
DCPA01-78-C-0279,  Work  Unit  2563F  (September  1979). 

6.  Carrier,  G.,  Fendell,  F.,  Feldman  P.,  and  Fink,  S.,  "Forced  Convection 
Extinction  of  a  Diffusion  Flame  Sustained  by  a  Charring  Body,"  Final  Report, 
Contract  EMW-C-0371,  FEMA,  WU  2563G,  TRW  Defense  and  Space  Systems  Group, 
Redondo  Beach,  CA  (May  1981). 

7.  Kanury,  A.M.,  Introduction  to  Combustion  Phenomena,  Gordon  and  Breach  Science 
Publishers,  New  York,  NY  (1975). 

8.  Williams,  F.A.,  Combustion  Theory,  Addison  Wesley,  Reading,  Mass.  (1965). 

9.  Bhagat,  P.M.,  "Wood  Charcoal  Combustion  and  the  Effects  of  Water  Application," 
pp.  275-291,  Combustion  and  Flame,  37,  (1980). 

10.  Kuchta,  J.M.,  Kant,  A.,  and  Damon,  G.H.,  Industrial  and  Engineering  Chemistry, 
44,  pp.  1559-1563,  (1952)  as  discussed  in  Appendix  P  of  Combustion  of 
Pulverized  Coal  by  M.A.  Field,  et  al ,  Published  by  the  British  Coal 
Utilisation  Research  Board,  Leatherhead,  England,  (1967). 

11.  Smith,  D.F.,  and  Gudmundsen,  A.,  Industrial  and  Engineering  Chemistry,  23, 
pp.  277-285,  (1934)  as  discussed  in  Appendix  P  of  Combustion  of  Pulverized 
Coal . 

12.  Parker,  A.S.,  and  Hottell ,  H.C.,  "Combustion  Rate  of  Carbon,"  pp.  1334-1341, 
Industrial  and  Engineering  Chemistry,  28(11) ,  (November  1936). 

13.  Martin,  S.B.,  Personal  Communication,  (1982) 


43 

14.  Kanury,  A.M.,  "Rate  of  charring  Combustion  in  a  Fire,"  pp.  1131-1142,  Four¬ 
teenth  Symposium  (International)  on  Combustion,  The  Combustion  Institute, 
Pittsburgh,  PA  (1973). 

15.  Kanury,  A.M.,  "Thermal  Decomposition  Kinetics  of  Wood  Pyrolysis,"  pp.  75-83, 
Combustion  and  Flame,  18,  (1972). 


rrvTv^ 


■J  jr  7”  ■  7*  t  / 


A  % 


O 

5  O 

s « 

111  u 
|T  III 
141 
Ul  U. 

—  141 

H  > 

°  P 

uj  *Z 

>  2 

z  a 

o  < 

o  z 


III 

14. 

• 

z 

< 


to 

m 

< 

Z 


tr 

UJ 

u. 

« 


< 

>  « 
O  « 

z  o 


60 
69 
O 
-1 

r  w  ^  . 

r  —  uj  r 


69 

CO 

< 

z 

z 

UJ 

0 

> 

X 

o 


2 

UJ 

5 

UJ 

0. 

60 

H 

O 

D 

8 


UJ 

z 


Ul 


UJ  X 
—  UJ 

•t  > 
O  r 
uj  J 

>  5 
z  c 
O  < 
o  Z 


©Q  ©@  ®@  © 


o 

a 


3 


ialfa 


/’  •  *  .  *  »'*•  , 


."•  • 

te.  ’  -  ' 


Figure  1.  A  Pool  rire  (a)  in  Huiescent  Air  and  (b)  in  a  Wind 


figure  4.  Peaction  Temperature  Dependence  on  Damkohler  Number. 
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Peaction  Temperature  Pependenco  on  Pamkohler  Number  When  ''ctivation  Energy  is  Low 


Reaction  Temperature  Dependence  on  Damkohler  Number  When  Activation  rneroy  is  Hioh. 


Map  of  Limiting  Damkohlcr  Number  for  Ignition  and  Extincti 
conditions:  Inverse  o*  Blast  Strength,  P  j,e  versus  Fire 
Strength,  q  . 


Figure  11.  Mechanism  of  Extinction  and  Peignition  of  Flanes  on  Charring  Solids. 
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Figure  13:  The  SRI  Shock  tube  Test  Section  Detail 


sticks  3/4-inch  square 
spaces  3/4-inch 


Designs:  (a)  3/4-inch  Stick  Crib,  and  (b)  Detail  of  3/8-inch  Stick  Crib. 


Figure  15:  Extinction  (  filled  symbols  )  -  No  Extinction  (  open  symbols  ) 
Coorrelation  for  Class  B  Fires  with  No  Barriers. 
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rinure  16.  Test  of  the  Hypothesis  for  the  E"ect  of  Barricade  on  Mess 


SUHHARY  OF  HEXANE  TEST  RESULTS 


**  These  tests  are  rereported  in  Ref.  2.  See  Table 


(Table  B.l  of  Ref.  2) 
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tains  the  blast.  The  fire  strength  is  described  for  liquid  sustains  the  blast.  The  fire  strength  is  described  for 
fuels  primarily  by  the  heat  of  combustion;  and  for  wood,  liquid  fuels  primarily  by  the  heat  of  combustion;  and  for 

it  is  mainly  described  by  the  preburn  time.  The  concept  wood,  it  is  mainly  described  by  the  preburn  time.  The 
is  substantiated  Dy  the  SRI  shock  tube  data  concept  is  substantiated  by  the  SRI  shock  tube  data. 
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